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FOREWARD 


This  document  is  a  supplement  to  the  Design  Notebook  for 
Naval  Air  Defense  Simulation*  SSP-G00-06-0049-82  dated  15  September 
1982fVThe  supplement  Is  a  detailed  description  of  the  AAW  Surveillance 
Simulation  which  was  originally  developed  for  OP-654  and  was  expanded 
under  the  current  contract  to  model  surveillance  by  radar  as  well  as 
ELINT  satellites  of  naval  and  air  forces.  Base  watching ,  various  forms 
of  barriers  and  flight  following  tactics  have  been  provided.  Space 
based  surveillance  vehicles  are  limited  to  circular  orbits. _ 
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INTRODUCTION 


TRWs  Waterwheel  Program  office  has  been  developing  the  COVERAGE 
simulation  for  OP-654,  OP-96,  other  customers,  and  under  company  IR&D  funding 
for  several  years.  Its  current  purpose  and  capability  Is  to  simulate 
detections  of  aircraft,  by  active  and  passive  satellites  and  OTH  radars.  The 
COVERAGE  simulation  requires  several  additional  programs  as  pre-processors  and 
post-processors.  The  purpose  of  this  documentation  Is  to  describe  Its 
capabilities  and  purposes.  Its  relationship  to  other  programs,  and  to  give  a 
sufficiently  detailed  description  of  the  code  to  permit  computer  programmers 
to  modify  the  code. 

At  this  writing,  the  COVERAGE  simulation  Is  undergoing  modifications 
as  a  part  of  current  Waterwheel  studies.  In  some  cases,  especially  for 
subroutines  BRBLD  and  BRRIER,  the  documentation  describes  modifications  about 
to  be  made  to  the  code  rather  than  the  current  code  which  will  shortly  be 
obsolete. 
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1. 


OVERVIEW 


1.1  BRIEF  DESCRIPTION  OF  COVERAGE  SIMULATION 

(U)  COVERAGE  Is  a  time-stepped  Monte  Carlo  simulation  which  models 
detections  of  aircraft  and  ships  by  active  and  passive  satellites.  A  high 
level  flow  diagram  Is  shown  In  Figure  1-1.  Satellite  and  target  (i.e., 
aircraft  and  ship)  motion  Is  modeled  explicitly.  At  each  time  step  the 
position  of  the  satellites  and  targets  are  updated  and  tests  are  performed 
to  see  If  any  of  the  satellites  detect  any  of  the  targets.  The  tests  for 
passive  detection  are: 

o  Is  the  target  emitting? 

o  does  the  emitter  Illuminate  any  spacecraft? 

o  Is  a  ground  station  visible  to  the  Illuminated  spacecraft? 

The  tests  for  active  detection  are: 

o  Is  the  target  emitting? 

o  does  the  emitter  Illuminate  any  spacecraft? 

o  Is  a  ground  station  visible  to  the  illuminated  spacecraft? 

o  Is  the  target  In  the  spacecraft's  detection  annulus  (of  which 
the  outer  edge  Is  power  limited  and  the  Inner  edge  Is  clutter 
limited)? 

o  for  aircraft  targets,  does  the  doppler  return  from  the  aircraft 
exceed  the  return  from  the  background  clutter? 

o  Is  the  radar  beam  pointed  at  the  target?  (Note  that  this 
simulation  models  mechanical  scanning  only  and  not  phased 
array.) 

o  a  random  draw  against  the  probability  of  detection. 

Active  and  passive  satellites  can  be  modeled  simultaneously  so  that,  for 
example,  an  active  satellite  can  be  cued  to  look  at  a  target  that  was 
Initially  detected  by  a  passive  satellite. 

(U)  In  each  replication  of  a  Monte  Carlo  run,  the  targets  commence 
their  courses  at  a  random  time,  so  that  the  satellites  are  In  a  different 
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position  when  the  bombers  take  off  and  when  ships  arrive  In  an  area  of 
Interest. 

* 

(U)  The  simulation  requires  three  Input  files  and  generates  three 
types  of  output  (see  Sections  1.3  and  1.4).  The  Input  files  describe: 

o  satellite  and  target  characteristics  (COVIN) 

o  target  courses  (CRS) 

o  density  of  background  aircraft  traffic  (BGTFIL) 

The  output  files  describe: 

o  time  of  detection  and  target  position  when  detected  (OBS, 
DETIM,  SCAN,  COVECH) 

o  satellite  hardware  performance  (SCAN  and  COVHST) 

o  data  for  input  to  a  graphics  program  (SCPLT,  TGPLT,  OVPLT). 

(U)  COVERAGE  is  used  In  conjunction  with  other  programs  (see 
Figure  1-2).  They  are: 


PROFIL  - 
CRSGEN  - 
STATCV  - 


WRNCOV  - 


AAWSIM  - 


MAPPRO  - 


generates  CRS  Input  files  for  COVERAGE  and  AAWSIM 

creates  a  graphics  input  file  from  a  CRS  file 

processes  DETIM  to  generate  statistics  on  coverage 
gaps 

processes  OBS  to  compute  probability  of  Initial 
detection  as  a  function  of  target  distance  along 
course 

a  one-on-one  Monte  Carlo  simulation  using  CRS  and 
OBS  to  model  engagements  between  an  Interceptor 
assisted  by  early  warning  aircraft  and  a  bomber  in 
which  the  Interceptor  and  early  warning  aircraft 
are  launched  and  vectored  based  on  external 
surveillance  reports. 

a  graphics  program  using  DISSPLA  subroutines  used 
to  draw  snapshots  of  COVERAGE  and  AAWSIM 
scenarios. 


(U)  All  these  programs  are  written  In  FORTRAN.  The  most  current 
versions  of  all  these  programs  compile,  load,  and  execute  on  CDC  computers 
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Figure  1-2.  (U)  Relationship  of  Coverage  to  Other  Programs  (U) 


and  on  DEC  VAX/780  VMS  computers.  Earlier  versions  of  COVERAGE,  AAWSIM, 
and  COVSTAT  have  compiled,  loaded,  and  executed  on  a  MODCOMP  Classic  III 
computer. 


'  (U)  These  programs  have  been  under  development  for  several  years 

5 

with  funding  from: 

J 

|  o  the  Systems  Analysis  Division  of  the  Chief  of  Naval  Operations 

•  Program  Planning  Office  (OP-96) 

o  the  Strategic  and  Theatre  Nuclear  Warfare  Division  of  the  Chief 
:  of  Naval  Operations  Plans,  Policy  and  Operations  Office  (OP-65) 

j  o  contractor  IR&D  funds 

V 

;  1.2  PURPOSE  OF  THE  COVERAGE  SIMULATION 

'i 

(U)  The  overall  purpose  of  the  COVERAGE  simulation  and  related 
!  programs  is  to  assess  the  utility  of  candidate  surveillance  architectures 

f  to  Naval  and  Air  Force  anti-air  warfare,  primarily  against  the  Soviet 

Backfire  bomber,  and  to  naval  anti-surface  warfare,  primarily  to  support 
|  over-the-horizon  tarqeting  of  the  Tomahawk  anti-ship  missile.  Production 

!  runs  can  be  made  to  study  mission  utility  as  a  function  of: 


i 


o 

o 

0 

0 

0 

0 

0 

0 


the  number  of  satellites 
satellite  orbital  parameters 
satellite  performance  characteristics 
number  and  location  of  ground  stations 
C3  time  late 

tradeoffs  between  satellites  and  non-space  surveillance  systems 
barrier  placement  (active  satellites  only) 
tasking  (active  satellites  only) 


(U)  Passive  satellite  performance  Is  modeled  in  terms  of  the 
spacecraft -target  geometry  required  for  detection.  Active  satellite 
performance  is  modeled  in  terms  of  slew  angles  and  antenna  feed  motion 
required  for  mechanical  scanning,  energy  balance,  and  the  number  of 
targets  that  can  be  tracked  simultaneously. 


(U)  Non-space  surveillance  systems  that  can  be  modeled  (and  the 
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program  that  models  them)  are: 


o  Early  warning  aircraft,  e.g.,  the  Navy  E-2C  and  the  Air  Force 
E-3A  (AAWSIM) 

o  Interceptors,  e.g.,  F-14  and  F-15  (AAWSIM) 

o  High  frequency  Over-the-Horizon  Backscatter  (OTH  B)  radar  (to 
be  added  to  COVERAGE) 

o  DEW  Line  radars  (to  be  added  to  COVERAGE). 

(U)  Tasking  of  active  satellites  addresses  how  many  carrier  battle 
groups  and  how  many  bomber  approaches  can  be  adequately  covered 
simultaneously. 

(U)  Mission  utility  for  anti-air  warfare,  or  more  specifically. 
Naval  Anti-air  Warfare  (AAW)  and  Air  Force  Airborne  Tactical  Warning  (ATW) 
to  defend  North  America,  is  quantified  in  terms  of 

o  the  probability  of  initial  detection  of  a  raid  as  a  function  of 
distance  from  the  carrier  battle  group  or  from  CONUS  targets 
(quantified  using  COVWRN) 

o  the  probability  that  an  interceptor  can  visually  identify 

and/or  shoot  at  a  bomber  beyond  some  keepout  range  (quantified 
using  AAWSIM) 

o  coverage  gap  statistics  (quantified  using  COVSTAT) 

Mission  Utility  for  Naval  Anti -Surface  Warfare  can  be  quantified, 
currently,  in  terms  of  coverage  gap  statistics  and  probability  of  initial 
detection. 
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v**»  V  T-T 


^  .i.i.  ♦-«  v.-.t  ;wg  v,  V-V. 


CASE  CPU  TIME  *  160.3 
TOTAL  CPU  TIME  «  161.9 
FORTRAN  STOP 

‘SDEBUG- 1  -EX I TSTATUS »  is  ' JsSYSTEM-S-NORMAL *  no*m*l  successful  completion 
DBG>EXIT 


1.3 


INPUT 


1.3.1  BGT  FIL  Input  File 

The  BGT  FIL  Is  a  file  describing  worldwide  background  aircraft  traffic 
densities.  There  are  24  records,  each  record  giving  the  densities  at  000  HRS 

GMT,  0100  HRS  GMT,  . 2300  HRS  GMT  for  10°  x  10°  latitude-longitude  bins. 

Densities  are  In  number  of  aircraft  per  million  square  miles.  Each  record  has 
17  rows  corresponding  to  10°  latitude  Increments  from  ±  90°  to  -  90°  and  36 
columns  corresponding  to  10*  longitude  bins  from  0°  to  360° • 

OBS  Output  File 

The  observation  output  file  Is  written  on  (unit  12)  every  time  a 
sensor  detects  an  aircraft  target.  Data  written  are: 

Column  1  Time  of  detection  (hours)  since  simulation  start; 

Column  2  Number  of  sensor  making  the  detection; 

Column  3  Course  number; 

Column  4  Attack  (replication)  number; 

Column  5  Target  number  (row  on  TARGTS  array); 

Column  6  Target  latitude  (degrees); 

Column  7  Target  longitude  (degrees); 

Column  8  Time  (hours)  since  target  began  course; 

Column  9  Bearing  in  degrees  (relative  to  north)  from  target  to 

sensor; 

Column  10  Distance  (n.ml.)  of  target  from  destination; 

Column  11  1  -  target  Is  classified;  0  -  target  Is  not 

classified; 


1.3.2  COVIN  Input  File 


NATTAK 

Integer 

NCRS 

Integer 

m 

LSEED 

Integer 

Odd  number 

BRTST 

Logical 

Always  FALSE 

SIMARR 

Logical 

TRUE 

FALSE 

TLIM1 

Real 

Hours 

TLIM2 

Real 

Hours 

Rl,  R2 

Real 

Hours 

THRSH  Real  Knots 

BELEV  Real  Degrees 

BW  Real  Degrees 

PD  Real  0-10 

PC  Real  0-10 

TPRINT  Real  Hours 

TPRINT2  Real  Hours 


Number  of  replications; 

Number  of  aircraft  and  ship 
courses  to  be  read  In  from  CRS 
file; 

Random  number  seed; 

Obsolete; 

All  bombers  arrive  at  missile 
launch  points  simultaneously; 

All  bombers  take  off 
simultaneously; 

(Passive  sensor  only)  Emission 
duration  immediately  after 
aircrafts  take  off; 

(Passive  sensor  only)  Aircraft 
emission  duration  Immediately 
prior  to  missile  launch; 

(Passive  sensor  only)  Aircraft 
emitter  usage  ratio  between 
emissions  at  beginning  and  end  of 
flight.  Emitters  are  used  R1 
hours  out  of  every  R2  hours; 
(Active  sensor  only)  Minimum 
dlscernable  velocity; 

(Active  sensor  only)  Outer 
annulus  edge  at  which  aircraft 
barriers  are  searched  (earth 
grazing  angle); 

(Active  sensor  only)  Beamwldth; 
(Active  sensor  only)  Probability 
of  detection  for  single  dwell; 
(Active  sensor  only)  Probability 
of  classification  -  obsolete; 
Obsolete; 

Obsolete; 


IFL61 

Integer 

0 

Do  not  compute  rotation  of  node; 

¥  0 

Compute  rotation  of  node; 

I  HE  AO 

Integer 

0 

Suppresses  Input  echo  on  unit  6; 

1 

Prints  input  echo  on  unit  6; 

I  PLOT 

Integer 

0 

Suppresses  output  on  SCPLT  and 
TGPLT  files  (in  main  program); 

1 

Prints  output  on  SCPLT  and  TGPLT 
files  (in  main  program); 

IPLOT2 

Integer 

0 

Suppresses  output  on  DETIM  file 
(In  SETL) 

1 

Prints  output  on  DETIM  file  (in 
SETL); 

LAA 

Logical 

TRUE 

(Active  sensor  only)  Debug  output 
on  SRT,  AA,  and  CLM  arrays; 

FALSE 

Suppresses  debug  output; 

NSCAN 

Integer 

• 

Replication  for  which  SCAN,  CVPLT, 
SCPLT,  and  TGPLT  output  Is 
desired.  Zero  If  no  output  Is 
desired  for  these  files; 

TPRINT 

Real 

- 

Obsolete 

NSW 

Integer 

1 

Satellite  coverage  area  Is  a 
circle  (one  edge); 

2 

Satellite  coverage  area  Is  an 
annulus  (two  edges); 

ISW 

Integer 

1 

Method  of  Inputlng  satellite 
coverage  area  edges; 

Earth  central  angle  for  each  edge; 

2 

Earth  grazing  angle  for  each  edge; 

3 

Spacecraft  look  angle  (from  nadir) 
for  each  edge; 

TMAX 

Real 

Hours 

Time  for  each  replication  or 
attack;  should  be  longer  than  the 
longest  course  time; 

NSC1 

Integer 

1-8 

Number  of  type  1  sensors; 

LTYP1 

Logical 

TRUE 

Type  1  sensors  are  active; 

FALSE 

Type  1  sensors  are  passive; 

v-  »r»r-v*Tr-. 


3  DTIMl 

jj 

j  OESC 

Real 

Hours 

Time  Increment  for  Type  1 
spacecraft; 

Orbital  parameters  for  spacecraft 
or  dimensions  of  coverage  area  for 
OTH  radars  and  DEW  Line; 

Orbital  parameters  (Initial 
conditions); 

R  OESC(l,  J) 

Real 

Feet 

Length  of  semi -ml nor  axis; 

|  OESC(2,  J) 

Real 

0.-1.0 

Eccentricity; 

|  OESC (3,  J) 

Real 

Degrees 

Inclination; 

I  OESC{4,  J) 

Real 

Degrees 

Longitude  of  ascending  node; 

I  OESC(5,  J) 

Real 

Degrees 

Argument  of  perigee; 

|  0ESC(6,  J) 

Real 

Degrees 

Mean  anomaly; 

OTH  B/DEW  Line  Dimensions; 

|  OESC(l,  J) 

Real 

Degrees 

Latitude  of  coverage  center; 

|  0ESC(2,  J) 

Real 

Degrees 

Longitude  of  coverage  center; 

p  OESC (3,  J) 

Real 

Oegrees 

Heading  (from  north)  from  center 
of  coverage  area  axis  of  symmetry 
(or  bisector); 

1  OESC (4,  J) 

Real 

Degrees 

Half  angle  of  width  of  coverage 

area; 

|  OESC(5,  J) 

Real 

n.ml. 

Radius  to  Inner  edge; 

I  0ESC(6,  J) 

Real 

n.ml. 

Radius  to  outer  edge; 

1  ORBPAR(NSC, 

I  1-2.  ISW) 

Real 

Degrees 

Satellite  coverage  area  dimensions 
(See  ISW); 

1  BATCAP 

Real 

Kilowatt  Hrs 
(KWH) 

(active  satellite  only)  Maximum 
storage  capacity  of  battery; 

I  PS0LAR 

Real 

KW 

(active  satellite  only)  Average 
^  power  Input  due  to  Insolation; 

1  PSPC 

Real 

KW 

(active  satellite  only)  Average 
power  required  for  base  load  of 
spacecraft; 

1  PPLON 

Real 

KW 

(active  satellite  only)  Power 
required  with  payload  on; 

|  PPLOFF 

Real 

KW 

1-16 

(Active  satellite  only)  Power 
requled  with  payload  off; 

ssssssi 
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CLSWID 

Real 

Degrees 

(active  sensor  only)  Maximum  feed 
angle  from  boresight.  (rotational 
degrees  about  suborbital  point); 

TGTTRN 

Real 

Degrees 

Maximum  expected  aircraft  turn; 

TGTSPD 

Real 

Knots 

Maximum  expected  target  speed 
change; 

DCPKPW 

Obsolete; 

DUTCYL 

Obsolete; 

SRCHRT ( 3) 

Real 

Degrees /Minute 

(active  sensor  only)  Search  rate 

for  barrier  search  for: 

1  Aircraft  barrier  at  outer  edge 
of  annulus; 

2  Aircraft  barrier  at  inner 
edge  of  annulus; 

3  Ship  barrier  at  outer  edge 
only; 


TTURN 

Real 

Minutes 

(active  sensor  only)  Feed  turn 
time  during  search; 

TRKCYL 

Real 

Hours 

(active  sensor  only)  Track  cycle 
time  (minimizes  energy  usage); 

EFF 

Obsolete; 

REQPEN 

Real 

Ratio 

Feed  turn  penalty  for 
reacquisition; 

TIMTRK 

Obsolete; 

PURTIM 

Obsolete; 

NSCZ 

Integer 

2-20 

Number  of  type  2  spacecraft; 

LTYPZ 

Logical 

FALSE 

Type  2  spacecraft  must  be  passive; 

DTIMZ 

Real 

Hrs 

Time  step  for  type  2  spacecraft; 

NSMCON 

Integer 

- 

Number  of  ground  station  sets; 

ISMTAB 

Ground  station  set  specifications; 

ISMTAB(1,  J) 

Integer 

1-15 

Set  Number; 

ISMTAB(2,  J) 

Integer 

- 

Target  one  start  number; 

ISMTA8(3,  J) 

Integer 

- 

Target  one  end  number; 

ISMTAB(4,  J) 

Integer 

- 

Target  two  start  number 

(first  G/S); 
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ISMTAB (5,  J)  Integer 
ISMTAB(6,  0)  Integer 
ISMTA8(7,  J)  Integer 


Target  two  end  number 
(last  G/S); 

Spacecraft  start  number  (0  =  any 
S/C); 

Spacecraft  end  number  (0  -  any 
S/C); 


Notes  on  ISMTAB; 

1.  Target  numbers  refer  to  rows  on  TARGTS  array; 

2.  Ground  stations  should  be  entered  sequent! lally  at  the  end 
of  the  TARGTS  array; 

3.  Entries  1-3  are  used  for  passive  sensors  only  and; 

4.  If  active  satellites  are  modelled,  ISMTAB  refers  to  active 
satellite  ground  stations  only; 


I 


ITSTOR 

Integer 

NUMCT 

Integer 

NTRGTS 

Integer 

TARGTS 

Real 

Zero 

Always  1 


Never  use; 

Number  of  simultaneous  observation 
by  spacecraft  required  for 
detection; 

Number  of  targets  to  be  Input  to 
the  TARGTS  array; 


k»l 


Input  area  for: 
o  Aircraft 
o  Ships 

o  Aircraft  barrier; 
o  Ship  barrier; 
o  Ground  stations; 


TARGTS  Input  for  aircraft  and  ships; 

TARGT$(J,  1)  1  Aircraft  raid  count; 

1  Ships  raid  count; 

TARGTS (J,  2)  Course  number; 

TARGTS(J,  3)  (Passive  sensor  only)  Emitter  scan  width  In  azimuth; 

TARGTS(J,  4)  (Passive  sensor  only)  Emitter  elevation  limits  1  or  2; 
TARGTS(J,  5)  (Passive  sensor  only)  Emitter's  lower  elevation  limit  (grazing 
angle  degrees  above  local  horizon); 


TARGTS(<Jt  6)  (Passive  sensor  only)  Emitter's  upper  elevation  limit  (grazing 
angle  degrees  above  local  horizon); 

TARGTS(J,  7)  Course  start  time  (hours)  -  computed  by  program.  Input  a  zero; 
TARGTS(J,  8)  For  passive  sensors.  Input  time  (hours)  after  takeoff  at  which 
emitter  first  turns  on; 

For  active  sensors,  enter  any  negative  number;  it  will  be 
replaced  with  the  time  of  initial  detection; 

TARGT$(J,  9)  Time  (in  hours)  of  course  duration;  is  computed  by  program; 
input  a  zero; 

Aircraft  and  Ship  Barriers; 

Note:  Ship  barriers  must  be  circles  centered  at  a  carrier  battle  group; 
aircraft  barriers  can  be  circles  or  annular  sections  centered  anywhere; 


TARGTS( J,  1) 

Must  be  zero; 

TARGTS(J,  2) 

Latitude  (degrees)  of  barrier  center; 

TARGTS(J,  3) 

Longitude  (degrees)  of  barrier  center; 

TARGTS( J,  4) 

Inner  radius  (n.mi.)  of  barrier  (must  be  zero  for 
circular  barrier); 

TARGTS(J,  5) 

Outer  radius  (n.mi.)  of  barrier; 

TARGTS(J,  6) 

Bearing  (degrees  from  north)  of  annular  barrier 
axis  of  symmetry; 

TARGTS(0,  7) 

Barrier  ID.  Must  be  an  odd  number  from  1  to  19  for 
aircraft  barrier  and  any  integer  from  21  to  30  for 
ship  barriers; 

TARGTS( J,  8) 

Barrier  width  (degrees); 

TARGTS(J,  9) 

Barrier  type  (1  for  aircraft,  2  for  ships); 

Ground  Stations: 

TARGTS(J,  1) 

Must  be  zero; 

TARGTS(J,  2) 

Latitude  (degrees); 

TARGTS(J,  3) 

Longitude  (degrees); 

TARGTS(J,  4) 

Elevation  coverage  limits  (1  or  2); 

TARGTS(0,  5) 

Earth  Grazing  angle  (degrees)  of  lower  elevation 
limit); 

TARGTS(J,  6) 


Earth  Grazing  angle  (degrees)  of  upper  elevation 
limit; 

TARGTS(J,  7)  ID;  Must  be  any  negative  integer  (entered  as  a 

real); 

TARGT$(J,  8-9)  Not  used; 

Revision  data: 


I 

Integer 

<  NSC1 

Number  of  type  1  sensors  to  be 
modelled; 

J 

Integer 

£  NSC  2 

Number  of  type  2  spacecraft  to  be 
modelled; 

IDO 

Integer 

1  NSMCON 

Number  of  ground  station  sets  to 
be  modelled; 

K 

Integer 

<  NTARGTS  + 

NSMCON 

Number  targets  to  be  modelled; 

1.3.3  COURSE  FILE 

(U)  This  file  defines  all  bomber  flight  paths  to  be  simulated  in  terms 
of  state  vectors.  Each  state  vector  consists  of: 

1)  Time  of  course  change  (hours) 

2)  Lattitude  of  the  course  change  (degrees) 

3)  Longitude  of  the  course  change  (degrees) 

4)  Target  acceleration  after  course  change  (knots/hr.) 

5)  Target  velocity  after  course  change  (knots) 

6)  Time  since  last  course  change  (hours) 

7)  Distance  flown  since  last  course  change 
(nautical  miles) 

8)  Total  distance  flown  thus  far  (nautical  miles) 

9)  Target  heading  after  course  change  (degrees) 

,  10)  Distance  from  carrier  (nautical  miles) 

A  separate  program  (PROFILE)  generates  elements  one  through  nine  for  each  the 
state  vectors.  AAWSIM  computes  the  tenth  element  of  each  state  vector. 


1.4  OUTPUT 


1.4.1  DETIM  Output  File 


DETIM  ts  written  on  (unit  8)  by  subroutine  SETL  every  time  a  new 
contact  (or  detection)  begins  of  a  barrier,  aircraft,  or  (for  passive  sensors 
only)  a  ground  station.  Data  written  are: 


Column  1 
Column  2 
Column  3 


Column  4 


Column  5 


Target  number  (rows  or  TARGTS  array); 
Attack/replication  number; 

Time  (hours)  since  aircraft  began  course  or  since  (for 
ground  stations  and  barriers)  replication  began; 

Down  or  gap  time  (hours)  since  end  of  last  contact  and 
start  of  this  contact; 

0  -  target  has  not  been  classified; 

1  -  target  is  classified; 


viaaKSuelii 


"J 


1.4.2  TGPLT,  SCPLT,  and  CVPLT  Output  Files 

These  files  have  a  special  format  required  by  the  MAPPRO  graphics 
program.  This  format  is 

FILID 
data  rows 
0,0, 0,0 
FILID 


data  rows 
0,0, 0,0 


repeated  an  unlimited  number  of  times 


End -of -file 

The  format  for  each  data  row  Is  LONGITUDE,  LATITUDE,  symbol  number, 
heading  where  all  angles  are  In  degrees.  Symbol  numbers  are  the  following 
Integers: 

1  Causes  a  DISSPLA  symbol  to  be  drawn  for  each  data  row; 

0  A  line  connects  all  the  points  indicated  on  the  data 

rows ; 

-  1  Do  not  use; 

-  2  Do  not  use; 

-  3  Draws  Backfire  symbol; 

-  4  Draws  F-14  symbol; 

-  5  Draws  CV  symbol; 

-  6  Draws  OTH  B  symbol; 

-  6  Draws  other  symbols  satellites,  RPV,  E-2C,  or  an  "X"; 


a 


FILIDs  must  be  ascending  real  numbers.  Free  format  is  used. 

TGPLT  writes  a  data  row  to  show  the  current  position  of  every  bomber 
at  every  simulator  time  step.  The  FILID  Is  the  simulator  time  at  which  the 
data  row  Is  written. 

SCPLT  writes  a  data  row  at  every  time  step  for 

o  suborbital  point  of  each  satellite; 
o  boreslght  almpolnt  for  each  cluster  (active 
satellite  only) 


o  Intersections  of  annulus  edges  with  aircraft  and 
ship  barrier  edges; 

Again,  the  FIL1D  Is  the  simulator  time  at  each  time  step. 

CVPLT  writes  scenarlc  data  rows.  Data  Is  Indicated  by  the  FILID. 

FIL1D  DATA 

1-50  Aircraft  and  ship  courses; 

51  -  100  Destination  of  aircraft  and  ships; 

101  -  150  Barrier  centers; 

200  Aircraft  and  ship  positions  when  detected; 


1.4.3  SCAN  Output  Files 


The  SCAN  file  shows  the  status  of  each  sensor  at  each  time  step, 
especially  by  indicating  which  targets  are  observed.  Data  written  for  each 


sensor  are: 


Column  1  Time  (hours); 

Column  2  For  each  target  In  TARGTS  array  o  if  not  detected; 

9  If  aircraft  cannot  be  detected  because  velocity 
component  along  line  of  sight  is  less  than  the  minimum 
discemable  velocity  (active  satellite  only); 

7  aircraft  (ship)  not  detected  because  probability  of 
detection  is  less  than  one  (active  satellite  only); 

A-Z  indicates 

o  barrier  aircraft,  or  ship  is  detectable; 
o  cluster  ID; 

Column  3  #  Indicates  a  scheduling  conflict 

blank  Indicates  no  scheduling  conflict  (active 
satellite  only) ; 

Column  4  letter  indicates  which  cluster  was  scheduled; (The 

remaining  data  are  for  active  satellites  only.) 

Column  5  Boreslght  slew  angle  since  last  time  step  (rotation 

angle  degrees  about  boreslght); 

Column  6  Off-boreslght  feed  angle  required  to  scan  cluster 

(rotation  angle  degrees  about  suborbital  point); 

Column  7  Aircraft  track  count  (the  number  of  aircraft  being 

tracked  at  current  time  step); 

Column  8  Battery  state  charge  (0  to  1005t)  ; 


Also,  there  Is  a  summary  column,  which  shows  for  each  target. 

8  target  Is  detected 
0  target  Is  not  observable 

018  target  Is  observable  by  I  satellites  but  Is 
not  detected  because  of  scheduling  conflicts; 
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2.0  OVERVIEW  OF  MAIN  PROGRAM 

This  chapter  describes  the  functions  performed  by  the  main  program 
code.  The  code  Is  referred  to  In  terms  of  the  FORTRAN  DO/LOOP  and  other 
labels  and  subroutines.  The  main  program  Is  long  and  performs  much  of  the 
simulation  logic.  The  overall  structure  of  the  main  program  (called  COV) 
Is: 


Declaration  Statements 
Comments  describing  Inputs 
Data  statements 
FORMAT  statements 

Input  and  Initialization  (labels  10-200) 

Time  loop  (labels  200-1000) 

Sensor  loop  (labels  295-600) 

Target  detection  loops  (labels  380-500) 


Target  detection  loops  Include  logic  or  subroutine  which  test  for  the 
following  sensor  -  target  combi nations.  Under  the  LOCATION  heading,  names 
refer  to  subroutines  and  numbers  refer  to  FORTRAN  labels  and  DO  LOOPs  in' 
the  main  program. 


SENSOR 

active  satellites 


passive  satellites 


TARGET 

aircraft  barriers 
aircraft 

ship  barriers 
ships 

ground  stations 

aircraft,  ships, 
and  ground  stations 


LOCATION 

BARRIER  (in  LOOP  385) 
BARRIER,  LOOP  500,  and 
828-855 
BARR IRS 

BARRIRS,  LOOP  500,  and 
828-855) 

355-380 

LOOP  500 


OTH,  B  radar  & 
DEW/LINE 


aircraft 


FXSENS 


Other  Important  sections  of  code  In  COV,  and  the  associated  sensor.  If 
applicable,  are: 


SENSOR 

FUNCTION 

LOCATION 

Calculate  start  &  end 
time  for  each  replication 
and  takeoff  times  for 
aircraft 

200-210 

all  satelltes 

Solar  eclipsing 

220-230 

ECLPSE 

Update  position  of  ships 
and  aircraft  at  each  time 
step 

260-295 

all  satellites 

Update  spacecraft  position 
(In  terms  of  suborbital  point) 

300-320 

CORBIT 

active  satellites 

Test  for  ground  station 
visibility 

355-380 

active  satellites 

Identify  candidate  search 
areas  at  which  boreslght 
should  be  pointed 

FSTEP 

500-580,  Including 
CLUSTR 

passive  satellites 

Test  for  ground  station 
visibility 

LOOP  625 

active  satellites 

Decide  at  which  search  area 
to  point  the  boreslght  of 
each  spacecraft 

CLSTPP,  just 
after  675 

— 

Output  SCAN  and  plotting  files 

900-997 

The  following  test  describe  COV  sequentially  from  label  10  to  the  end  In 
greater,  but  not  complete,  detail: 


LABEL/LOOP/SUBROUTINE  DESCRIPTION 

LOOP  10  Rewind  files. 

COMBLK  Initialize  earth  gravitational  constant 

(GMU),  earth  rotation  rate  (ME)  and  earth 

radius  (RE)  In  distance,  time,  angle  units  of 
feet,  hours,  and  degrees. 


SECOND 


Save  current  value  of  CPU  time 


ZERRV 

ZERIV 

20-30 

40-45 

45-50 


LOOP  75 


LOOP  87 
LOOP  90 


LOOP  100 


105-140 


LOOP  155 


LOOP  165 


165-185 

BRBLD 

BRBLDS 

UPR1 

JDORGF 

LIHTRP 


Set  arrays  to  zero. 


Read  COVIN  file  (describes  sensors  and 
scenarios). 

Echo  COVIN  file  on  COVOUT. 

Determine  time  step  (DTIME)  and  IRATIO,  which 
makes  It  possible  toupdate  and  test  the  two 
types  of  satellites  at  different  time  steps 
(DTIM1  and  DTIM2).  TIME  Is  subsequently  set 
to  zero. 

Read  the  CRS  file  (describes  routes  followed 
by  aircraft  and  ships)  and  compute  distance 
between  starting  point  of  each  aircraft  or 
ship  course  and  the  location  of  the  carrier 
or  base  that  the  aircraft  or  ship  Is 
attacking. 

Write  aircraft/ship  courses  on  COVPLT  file. 

Write  locations  of  carriers  or  bases  that 
aircraft  or  ships  are  attacking  on  COVPLT 
file. 

Echo  CRS  file  on  COVSUM  file. 

Write  header  for  SCAN  file.  Contents  of 
header  Is  determined  by  the  number  of  sensors 
and  the  number  of  targets  (l.e.,  aircraft  and 
ship  barrier,  alrcrft,  ships,  and  ground 
stations. 

Spacecraft  Initialization.  DNODEH  Is 
procession  of  the  mode  In  degrees.  XM  Is 
spacecraft  velocity  and  HA  and  HP  are  apogee 
and  perigee. 

Converts  some  entries  of  the  TARGTS  array 
from  real  to  Integer. 

8RBLD  and  BRBLDS  set  up  arrays  that  are 
referred  to  In  processing  aircraft  and 
ship  barriers.  UPR1  a  uniform  random 
number  generator  and  Is  called  with 
JDORGF  to  compute  a  random  starting  date 
for  eclipsing.  Remaining  code  Initializes 
variables  for  computing  active  satellite 
energy  balance.  LINTRP  Is  a  linear 
Interpolation  routine. 


E«T>;,rv  --  .*-  «. 


Jwlvl 


2-3 


185-200 

200 

200  to  220 

220-230 

BGTINP 

ECLPSE 

235-250 

250-295 

260-263 

BGNCRS 

ENDCRS 

263-275 

275-295 

LOOP  500 


Compute  TMAX,  the  time  at  which  the 

simulation  ends. 

Start  of  time  loop,  in  which  sensor  and 
target  positions  are  updated  and  tests  are 
made  to  see  if  any  sensors  detect  any 
targets. 

For  each  attack  (l.e.,  replication),  compute 
simulated  time  at  which  this  replication 
begins  (SATTAK)  and  ends  (TATTAK)  and  at  what 
time  each  of  the  aircraft/ships  begins  its 
course  (TARGTS  (NTAR.7)). 

BGTINP  updates  clock  and  reads  in  new 
aircraft  background  density  data  If  hour 
has  changed.  ECLPSE  checks  for  solar 
eclipsing  of  satellite. 

Initialization  required  each  time  the  time 
loop  Is  updated. 

Update  position  of  moving  targets  and  DELTIM 
of  fixed  and  moving  targets.  DELTIM  Is  the 
difference  between  the  current  TIME  and  the 
time  the  replication  began  (for  fixed 
targets)  or  the  time  the  aircraft  or  ship 

began  Its  course  (for  moving  targets). 

Check  If  aircraft/ship  has  begun  or  ended 
course.  BGNCRS  Is  called  just  before  an 
aircraft  begins  Its  course  so  that  coverage 
statistics  will  Include  the  detection  or  gap 
occurring  when  the  aircraft  begins  Its 

course.  ENDCRS  Is  called  for  the  same 
purpose  just  after  an  aircraft  ends  Its 
course.  Otherwise,  If  for  example,  an 

aircraft  course  took  three  hours  and  the 

first  detection  did  not  occur  until  two  hours 
Into  the  flight,  that  two  hour  gap  would  not 
be  Inlucded  In  CVSTAT  statistics  If  BGNCRS 
were  not  used. 

Interpolation  between  two  state  vectors  from 
the  CRS  bit  to  get  aircraft/ship  position, 
speed,  and  heading. 

When  target  happens  within  six  minutes  from  a 
state  vector,  position,  course,  and  speed  are 
obtained  directly  from  the  state  vector 
Instead  of  from  Interpolation. 

Sensor  loop.  Performs  spacecraft,  OTH  B,  and 
DEW  Line  tests  for  current  time. 


2-4 


295-300 


300-320 

CORBIT 

320-350 


355-370 

FXSENS 

LOOP  385 
BRRIER 

LOOP 


LOOP  500 


405-425 

425-440 


Ships  sensor  loop  when  needed  if  the  two 
types  of  spacecraft  are  tested  at  different 
time  increments  (DTIM1  and  DTIM2). 

Update  spacecraft  suborbital  point  for 
current  time  step. 

Compute  earth  central  angle  for  inner  and 
outer  edge  of  spacecraft  detection  annulus 
based  on  spacecraft  altitude,  and  satellite 
look  angle  from  nadir  (BETA)  or  earth 
elevation  angle  above  horizon  (EPS). 

Test  if  active  satellites  have  a  ground 
station  visible.  If  not,  skip  rest  of  sensor 
loop. 

If  current  sensor  Is  an  OTH  B  radar  or  the 
DEWLSC  line  barrier,  test  for  detection  of 
aircraft  and  then  ship  to  end  of  sensor  loop. 

Test  ifany  aircraft  barriers  can  be  searched 
at  current  time  by  current  spacecraft 
active.  If  it  Is,  store  boundaries  of  search 
area  in  SRT  array  for  spacecraft  scheduling. 

Test  if  any  ship  barriers  can  be  searched  at 
current  time  by  current  spacecraft,  if 

ctlve.  If  so,  store  search  area  boundaries 
on  SRT  array  for  spacecraft  scheduling. 

For  active  satellites,  tests  for  detection  of 
aircraft  that  have  already  been  initially 
detected,  either  in  a  barrier,  by  a  passive 
satellite,  by  the  DEW  line,  or  by  an  OTH  B 

radar.  For  passive  satellites,  test  for 

detection  of  a  ship  or  aircraft  and  for 

ground  station  visibility. 

Retrieve  target  position,  course,  speed  and 
DECTIM  from  TPOS. 

Test  Is  target  Is  in  annulas  by  comparing 
earth  central  angles  between  suborbital  point 
and  target  and  suborbital  point  and  nular 
edges. 


2-5 


445-465 

SETL 


465-475 

JMEQ 

FSTEP 


500 

500-565 


Complete  detection  tests  for  passive 
satellites,  namely:  is  aircraft/ship 

emitting?  If  emitting,  is  it  emitting  in  the 
direction  of  the  aircraft? 

Emitter  on/off  usage  is  described  by  TUMI, 
TLIM2,  Rl,  and  R2  (see  Table  ).  If  tests 
are  passed  and  ground  station  visibility  Is 
not  being  tested,  the  detection  is  written  as 
CVSCAN,  CVOBS,  and  SETL  is  called.  SETL 
updates  coverage  statistics  array  for  current 
replication  and  may  cause  write  on  CVTAP8 
(Input  to  CVSTAT).  If  ground  station 

visibility  is  being  tested, then  target 
observability  is  Indicated  on  the  IRFTAB 
array.  If  current  target  is  a  ground 
station,  then  its  observability  is  noted  in 
ISCAN  and  IRFTAB  and  SETL  is  called. 

Perform  doppler  test  for  active  detection 
of  aircraft  by  calling  JMEQ.  JMEQ  computer 
DDOP,  the  component  of  aircraft  velocity 
along  the  llne-of-sight  between  the  target 
and  the  spacecraft.  For  detection  to  occur, 
the  absolute  value  of  DDOP  must  exceed  the 
doppler  return  from  clutter,  which  Is  Input 
as  THRESH.  If  the  aircraft  Is  detectable, 
then  FSTEP  is  called  to  store  aircraft 
position  for  subsequent  use  in  deciding  where 
to  point  the  boresight  (In  subroutines  CLUSTR 
and  CLSTPP) . 

End  sensor  loop. 

By  now,  the  SRT  array  holds  the  spacecraft 
look  angle  (from  radar)  and  the  bearing 
(relative  to  spacecraft  velocity)  for  each 
target  and  for  each  search  area  end  point 
detectable  from  the  current  spacecraft.  The 
AA  array  is  now  computed.  It  stores  the 
suborbital-polnt-centered  angle  between  every 
pair  of  targets  (and  search  area  end  points) 
observable  by  the  current  spacecraft.  This 
Information  is  used  by  CLUSTR.  If  search 
area  endpoints  exist,  then  the  search  area  ID 
Is  stored  on  the  AA  diagonal,  which  otherwise 
would  be  zero. 


565-570 


Optional  listing  of  AA  for  debug  purposes. 


570-580 

CLUSTR 


580-630 

SETL 


650-675 


Computes  bores Ight  pointing  options  for 
current  spacecraft  by  partioning  detectable 
targets  and  search  area  into  clusters. 
Clusters  are  groups  of  targets  and  search 
areas  which  can  be  scanned  froma  fixed 
bores ight  position  by  moving  the  antenna 
feed.  After  CLUSTR  has  been  called  for  all 
spacecraft,  then  CLSTPP  is  called  (after 
label  678)  to  decide  where  to  aim  the 
boresight  of  each  spacecraft.  For  fixed 
sensors,  i.e.,  OTH  B  and  DEW  line,  all 
observable  targets  are  assumed  to  be  in  one 
cluster  and  the  CLUSTR  functions  are 
performed  in  FXSENS. 

Ground  station  logic  for  passive  spacecraft. 
The  ISMTAB  array  specifies  that  particular 
spacecraft  (ISMTAB(6-7,I))  can  report 
detections  of  particular  target 
(ISMTAB(2-3,I))  only  if  the  spacecraft  sees  a 
target  and  particular  ground  station 
locations  (ISMTAB(4-5,I))  simultaneously. 
For  targets  and  ground  stations,  ISMTAB 
refers  to  rows  on  the  TARGTS/ITRGTS  array. 
If  detection  occurs,  SETL  Is  called  to  update 
coverage  statistics  and,  perhaps,  write  on 
CVTAP8,  JSCAN  Is  set,  and  the  latitude  and 
longitude  of  the  observed  target  Is  written 
on  tape  8. 

ITSTOR  is  never  used. 


675-LOOP  815 
CLUSTR 

LOOP  815 
LOOP  875 


Decide  at  which  cluster  to  aim  active 
satellite  boresight.  This  process  Is  called 
scheduling.  CLSTPP  knows  the  clusters  of  all 
spacecraft. 

Update  'he  ISCAN  summary  (last  column)  to 
show  how  many  sensors  detect  each  target. 

Each  sensor  has  now  been  assigned,  or 
scheduled,  to  look  at  a  particular  cluster. 
This  loop  updates  active  spacecraft  hardware 
requirement  parameters  (state  of  change,  feed 
angle,  lew  angle,  and  number  of  aircraft  to 
be  tracked)  and  determines,  based  on  a  random 
draw  against  an  input  probability  of 
detctlon,  whether  a  particular  ship  or 
aircraft  Is  detected.  If  detection  occurs, 
ISCAN  Is  set,  CVOBS  Is  written  on,  and  SETL 
is  called. 


2-7 


Wf 


00875  -  D0855 
STATUD 


00855  -  830 
UPRI  (TSS) 
RAN  (VAX) 


830-855 

SETL 


855-870 

STATUD 


870-875 

ECLPSF 

PWRBAL 

STATUD 

875-900 


SCSCH  (3, IDO)  stores  the  number  of  the 
scheduled  clusters  (in  terms  of  a  row  on  the 
CLM/CLM2/CLM3  arrays).  STATUD  is  called 
twice  to  update  histogram  statistics.  The 
first  call  updates  the  number  of  aircraft 
that  must  be  tracked,  which  is  the  sum  of  the 
background  and  threat  aircraft  in  search 
areas  and  In  (smaller  search)  areas  where 
threat  aircraft  are  being  reacquired  and 
tracked.  The  second  call  updates  the  ratio 
of  the  number  of  aircraft  that  must  be 
tracked  (track  requirements)  to  the  numbe*  of 
aircraft  that  there  is  time  to  track  (t|*ack 
capacity).  A  message  is  written  whenever  the 
ratio  is  greater  than  one. 

For  each  aircraft  or  ship  target,  a  uniform 
random  draw  is  made  to  test  against  an  input 
probability  of  detection  (PD).  If  the  draw 
is  greater  than  PP,  detection  does  not  occur 
and  ISCAN  Is  set  to  -7,  which  shows  up  in  the 
SCAN  printout  as  a  +7. 

If  an  aircraft  is  being  detected  (in  a 
search  area)  for  the  first  time, the  time  of 
the  initial  detection  Is  stored  In  the  eighth 
element  of  the  TARGTS/ITRGTS  row  for  that 
target.  For  detected  targets  ISCAN  Is  set  to 
the  number  of  the  scheduled  cluster,  SETL  is 
called,  and  CVOBS  is  written  on. 

Based  on  the  boreslght  aim  points  at  the 
current  and  previous  time  steps,  spacecraft 
slew  angle  (centered  at  the  suborbital  point) 
is  computed  and  stored  in  CLBG  (5,5,  NSC)  and 
STATUD  Is  called  to  update  slew  angle 
histogram  statistics.  STATUD  Is  also  called 
to  update  histogram  statistics  on  cluster 
width  (the  suborbital  point  centered  angle 
between  the  targets  or  search  area  corners  at 
the  edges  of  a  cluster).  j 

Update  power  balance  statistics  for  current 
spcecraft.  ECLPSF  is  called  to  determine  If 
current  spacecraft  Is  eclipsed.  In  which  case 
solar  panels  are  not  charged. 

For  targets  currently  being  detected,  set 
time  of'  last  detection ( TARGTS( NTAR ,  11 ) )  to 
current  time. 
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90-995 


Write  on  output  files  CVSCAN,  CVPLT,  CVTG, 
and  CVSC  If  DOUT  Is  true.  DOUT  prevents 
writes  at  the  smaller  time  step  in  DTIM1  and 
0TIM2  are  different. 


900-920 

CTOALF 


920-955 


955-985 


985-995 


995-1000 


1000-1050 

STATPT 

SUMMRY 


Write  on  CVSCAN  file.  CTOALF  converts 
positive  numbers  stored  in  ISCAN  to  alphabet 
letters.  Different  SCAN  formats  are  used 
depending  on  whether  one  or  two  types  of 
spacecraft  are  being  simulated.  CVSCAN 
shows,  for  each  sensor,  which  targets  are 
detected,  the  different  clusters,  which 
pluster  Is  scheduled,  and  hardware 
requirements  parameters. 

Writes  longitude  and  latitude  of  current 
position  of  detected  moving  targets  on  CVPLT 
file. 

Writes  longitude,  latitude,  and  heading  of 
suborbital  points  of  spacecraft  on  CVSC 
file.  If  active  spacecraft  are  being 
simulated,  the  longitude  and  latitude  of 
boreslght  aim  points  are  written.  If  active 
spacecraft  are  simulated  and  their  annulus 
edges  Intersect  barriers  (thus  creating 
search  areas),  the  longitudes  and  latitudes 
of  these  Intersections,  which  jre  the  search 
area  endpoints,  are  also  written  on  SCPLT. 

Write  longitude,  latitude,  and  bearing  of  all 
moving  targets  at  current  time  step  on  TGPLT 
file. 

End  of  time  loop.  If  TMAX  is  exceeded,  exit 
loop  and  end  simulation.  Otherwise, 
Increment  the  number  of  time  steps  (ITIME), 
time,  clock,  and  date. 

Call  final  statistics  subroutines,  write 
endflles,  rewind  files,  and  stop.  STATPT 
writes  histogram  dates  on  tape  4  and  SUMMRY 
computes  coverage  statistics  for  last 
replication  only  and  writes  then  on  CVECHO. 


whs 


3.  SUBROUTINES 

3.1  INPUT 


V*  V 


3.1.1  SUBROUTINE  BRBLD 
PURPOSE: 

Reads  Input  data  about  aircraft  barriers  from  targets  array  and  stores 
it  in  BRREF/IBRREF  for  processing  aircraft  barriers  in  BRRIER.  Barriers  are 
modeled  as  circles  or  annular  sections. 

INPUTS: 

NTRGTT 

TARGTS  (NTAR.l) 

TARGTS  (NTAR.2) 

TARGTS  (NTAR.3) 

TARGTS  (NTAR.4) 

TARGTS  (NTAR.5) 

TARGTS  (NTAR.6) 

TARGTS  (NTAR.7) 

TARGTS  (NTAR,8) 

TARGTS  (NTAR.9) 


Number  of  input  rows  to  TARGTS  array  inputs  for 
aircraft  barriers; 

Must  be  zero  if  row  contains  barrier  information; 
Longitude  of  barrier  center  (deg); 

Latitude  of  barrier  center  (deg); 

Radius  to  inner  edge  of  barrier  (n.mi.); 

Radius  to  outer  edge  of  barrier  (n.mi.); 

Bearing  of  bisector  of  barrier  (deg.  from  north); 
Fixed  target  ID;  for  aircraft  barriers  must  be  an 
odd  number  from  1  to  19; 

Barrier  width  (half  angle); 

Not  used 


For  circular  aircraft  barriers,  TARGTS  (NTAR,4)  must  be  zero  and  the  sixth  and 
eighth  entries  are  not  read.  Annular  barriers  must  have  two  rows  on  the 
TARGTS  array,  the  first  of  which  must  have  the  Inputs  described  above  and  the 
second  of  which  must  have  the  first  entry  set  to  zero  are  the  seventh  set  to 
the  same  barrier  ID. 


TARGTS  array  inputs  for  aircraft. 


TARGTS  (NTAR,1) 
TARGTS  (NTAR.2) 
TARGTS  (NTAR.3) 

TARGTS  (NTAR.4) 
TARGTS  (NTAR.5) 

TARGTS  (NTAR.6) 
TARGTS  (NTAR.7) 


Must  be  greater  than  zero; 

Course  number  In  CRS  input  file; 

Not  used 
Not  used 

Not  used 
Not  used 

Time  after  start  of  replication  that  bomber 
begins  course.  Is  not  an  Input  but  is  computed  by 
program; 
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TARGTS  (NTAR,8) 


Set  to  any  negative  number.  If  and  when  aircraft 
Is  detected  In  a  barrier,  this  entry  Is  set  to 
the  time  at  which  Initial  detection  occurs. 

TARGTS  (NRAT.9)  Length  of  flight  (hours); 

OUTPUTS 

NUMBR  Number  of  aircraft  barriers  in  TARGTS  array; 

IBRREF  (l.IBR)  Barrier  ID; 

IBRREF  (2.IBR)  Number  of  (first)  row  in  TARGTS  array  with 

Information  about  barrier  center,  radii,  and 
bearings; 

IBRREF  (3,I8R)  Number  of  (second)  row  in  TARGTS  array,  with 

barrier  ID  only; 

IBRREF  (4,IBR)  Latitude  of  pole  of  a  great  circle  going  through 

the  barrier  center  and  defining  one  side  of  the 
barrier; 

IBRREF  (5,IBR)  Longitude  of  pole  first  described; 

IBRREF  (6,IBR)  Latitude  of  pole  of  great  circle  of  radial 

defining  other  side  of  barrier; 

IBRREF  (7,  IBR)  Longitude  of  pole  just  described; 

Has  a  shorter  current  detection  period  than  that 
of  any  of  the  targets  seen  by  the  other  cluster; 
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METHOD: 


Aircraft  barriers  are  defined  as  annular  sectors: 


This  subroutine  searches  the  TARGTS  array  to  build  the  BRREF  array.  The  poles 
for  the  radlals  are  computed  by  adding  90  to  the  bearing  of  a  radial  and 
finding  a  point  1/4  of  an  earth  radius  away  from  the  center  along  that 
bearing.  If  the  inner  radius  is  zero,  then  a  circular  barrier  is  assumed  and 
no  poles  are  calculated.  The  data  described  here  are  used  In  subroutine 
BRRIER. 


3.1.2  SUBROUTINE  BRBLDS 


PURPOSE: 

Sets  up  a  reference  array,  BREFS,  for  use  by  subroutine  BRRIRS  In 
processing  ship  barriers.  Like  BRBLD,  which  Is  for  aircraft  barriers,  BRBLDS 
gets  the  data  for  the  reference  array  by  searching  the  TAROTS  array.  Ship 
barriers  are  models  as  circles  centered  at  a  carrier  battle  group  with  an 
Input  radius  specified  on  the  TARGTS  array.  They  are  searched  only  at  the 
outer  edge  of  the  annulus.  Ship  barriers  are  used  to  model  radar  satellite 
searches  of  enemy  ships  for  providing  over-the-horizon  targeting  of  Tomahawk. 

INPUTS: 


For  ship  barriers,  TARGTS  must  have  the  following  Inputs: 


TARGTS  (l.NTAR) 
TARGTS  (2,NTAR) 
TARGTS  (3.NTAR) 
TARGTS  (7.NTAR) 


Zero; 

Latitude  of  barrier  center; 

Longitude  of  barrier  center; 

Barrier  ID;  must  be  a  number  from  21  to  30 
(Implies  a  maximum  of  ten  ship  barriers); 


TARGTS  (8.NTAR)  Radius  of  barrier  (n.ml.); 


To  Input  ships,'  TARGTS  must  have: 


TARGTS  (NTAR.l) 
TARGTS  (NTAR.2) 
TARGTS  (NTAR.7) 

TARGTS  (NTAR.8) 

TARGTS  (NTAR.9) 


A  negative  number; 

Number  of  course  Input  In  CRS  file; 

Time,  (In  hrs)  at  which  ship  begins  Its  course 
(computed  by  program); 

Time  at  which  ship  Is  Initially  detected 
(computed  by  program); 

Duration  of  ship  course  (hrs); 


OUTPUTS: 

NUMBRS 

IBREFS  (1,IBR) 
IBREFS  (2, IBR) 
IBREFS  (3, IBR) 
IBREFS  (4, IBR) 


Number  of  ship  barriers  Input  In  TARGTS  array; 
Barrier  ID  number; 

Row  In  TARGTS  array  describing  this  barrier; 
Not  used 
Not  used 


IBREFS  (5,IBR) 
BREFS  (6,IBR) 
BREFS  (7.IBR) 
BREFS  (8,IBR) 
BREFS  (9,IBR) 


Not  used 

Latitude  of  barrier  center  (degrees); 
Longitude  of  barrier  center  (degrees); 
Barrier  radius  (n.ral.); 

Barrier  radius  (degrees); 


METHOD: 

As  part  of  program  Initialization,  BRBLDS  tests  every  row  of  the 
TAR6TS  array  for  a  zero  In  the  first  entry  and  a  number  between  21  and  30  In 
the  seventh  entry.  If  these  tests  succeed,  then  the  row  describes  a  ship 
barrier  and  the  appropriate  Information  Is  put  Into  the  BREFS  array. 


gjj 


•tv*: 
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3.2  SUBROUTINE  CORBIT 

Subroutine  CORBIT  calculates  the  suborbital  point  and  the  velocity  of 
the  suborbital  point. 

Spacecraft  position  and  velocity  in  geocentric  inertial  spherical 
coordinates  are  translated  Into  earth  geographical  coordinates.  The  elements 
used  to  describe  the  circular  orbit  for  an  earth  satellite  are  listed  in  table 
1. 

CORBIT  first  calculates  the  orientation  of  the  orbit  plane,  that  is, 
the  longitude  of  the  ascending  node  at  the  time  of  interest, 

ON  *  TIME(DNODE-WE)  +  00 

and  the  argument  of  perigee;  the  rotation  of  theorbit  within  the  orbital  plane 
PM  *  (XNU  *  TIME  +  0M)  MOO  360° 

where  MOO  Is  the  remainder  function. 

The  position  of  the  satellite  at  the  time  of  Interest  thus  calculated 
is  translated  to  earth  coordinates. 


Symbol 

TABLE  1 

Element 

Definition 

DI 

Inclination 

Angle  between  orbit 

XNU 

Mean  Angular  Motion 

plane  and  earth's 
equatorial  plane 

Angular  velocity  in  the 

A 

Distance 

orbital  plane  y/a3 
Distance  from  earth 

center  to  the  satellite 

DNODE 

Precession  of 

Angular  velocity  of  the 

ascending  node 

orbit  plane  about  the 

PM 

Argument  of  perigee 

line  normal  to  the 
earth's  equatorial  plane 
and  throuqh  the  origin 
(earth  center) 

The  rotation  of  the 

00 

Right  ascension  of 

orbit  within  the  orbital 
plane 

Longitude  of  ascending 

the  ascending  node 

node  at  reference  time 

OM 

Mean  anomaly  at  epoch 

360°  x  At/p  where  P  Is 

To 

Epoch 

the  period  and  t  Is  the 

time  difference  between 

the  epoch  and  the  last 
perigee  passage  before 
epoch 

Time  for  which  elements 

ON 

Longitude  of 

are  specified 

Orientation  of  the 

ascending  node 

orbital  plane  at  the 

WE 

Earth  rotation  rate 

time  of  Interest 

GMU 

(degrees/hour) 
Gravitational  constant 
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The  coordinate  of  the  point  is  (PM,  0)  in  Inertial  coordinates 

In  earth  reference  the  old  North  pole  is  at  (0N-90*,  90-DI),  the 
coordinates  of  the  point  in  the  new  system  0^=  90° 

COS (90  -  LAT)  =  COS(DI)  *  0  +  SIN(DI)COS(PM  -  90) 

SIN(LAT)  *  SIN(DI)  SIN  PM 
LAT  =  SIN’1  SIN( 01 ) SIN  PM) 

SIN (L0N-90)  =  SIN(90  -  PM)  SIN  0°/SIN(90  -  LAT) 

COS(LON)  *  COS ( PM ) /COS ( L AT ) 

LON  =  COS  ‘1  |C0S(PM)C0S(LAT)|  . 


3.3  ACTIVE  SENSOR  LOGIC 


3.3.1  SUBROUTINE  BGT 
PURPOSE: 

This  subroutine  is  responsible  for  computing  the  background  traffic 
in  a  specified  area. 

INPUT  ARGUMENTS: 

AREA,  ALAT,  ALON,  BGTDEN 

Three  calling  sequence  inputs  are  used  to  define  the  area  of 
interest,  and  a  fourth  variable,  the  background  density  array,  is  provided 
through  X8  Common. 

AREA  Area  of  region  in  square  miles; 

ALAT  Latitude  representative  of  region  in  degrees; 

ALON  Lonqitude  representative  of  region  in  degrees; 

BGTDEN  Background  traffic  density  data,  stored  in  an  18  x  36; 
array,  in  aircraft  per  million  square  miles. 

(BGTDEN  is  contained  in  COMMON  X8); 

OUTPUT  ARGUMENTS:  TRAFIC 

The  aircraft  count  resultina  from  the  computations  of  this  routine 
are  returned  via  the  calling  sequence  variable  TRAFIC. 

TRAFIC  Number  of  aircraft  in  the  region. 


METHOD: 

This  routine  computes  the  aircraft  count  in  a  region  based  upon  the 
area  and  the  density  associated  with  the  region.  The  region's  latitude 
and  longitude  are  converted  Into  array  subscripts  by  the  following 
formulae: 


* 


LAT  =  9  -  lFlX(l.E-6  +  ALAT/10. ) 

10  -  lFIX(l.E-6  +  ALAT/10.) 
LON  =  18  +  lFlX(l.E-6  +  AL0N/10.) 

19  +  1F1X( 1.2-6  +  AL0N/10.) 


if  ALAT  >  0. 
if  ALAT  <  0. 
if  ALON  1  0. 
if  ALON  <  0. 


The  traffic  count  Is  then  computed  by  the  following  statement,  which 
includes  the  conversion  from  aircraft  per  million  square  miles: 


V'J 


TRAFIC  *  BGTDEN(LAT.LON)  x  AREA/10. E6 
REFERENCES: 

The  routine  is  called  from  BGTSCN  and  CLUSTR. 


6 
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3.3.2  Subroutine  BGTINP 
PURPOSE: 

This  routine  is  responsible  for  keeping  the  background  traffic  data 
current  with  the  simulation  time.  The  background  traffic  file  B6TFIL 
contains  densities  by  hour  for  the  ten  degree  square  segments  covering  the 
surface  of  the  globe.  As  model  time  progresses,  the  density  data  must  be 
revised  for  each  hour. 

INPUT  ARGUMENTS: 

CLOCK,  BGTTIM,  BGTDEN 

The  routine  receives  the  simulated  clock  time  (24  hour,  with 
fractions  of  an  hour)  through  the  calling  sequence.  It  maintains  BGTTIM, 
the  time  associated  with  the  current  density  data  set,  and  BGTDEN,  the  18 
by  36  array  of  background  traffic  densities,  the  density  data  set.  These 
latter  two  items  are  input  from  the  BGTFIL  on  FORTRAN  logical  unit  number 
20  and  are  maintained  in  COMMON  X8. 

CLOCK  24.00  hour  simulation  clock  time 

BGTTIM  24.00  hour  time,  in  hourly  increments,  associated  with 
BGTDEN 

BGTDEN  Background  traffic  density,  stored  in  an  18  x  36  array, 
in  aircraft  per  million  square  miles. 

s  1 

OUTPUT  ARGUMENTS:  ' 

BGTTIM,  BGTDEN 

These  arguments  are  tn '  current  time  and  the  current  background 
traffic  density  data  as  updated  from  the  BGTFIL  as  necessary.  They  are 
defined  above. 


METHOD : 


The  CLOCK  time,  maintained  in  the  main  program,  is  compared  to  the 
BGTTIM  value.  Initially,  B6TTIM  is  set  to  -10.E20.  If  CLOCK  is  greater 
than  BGTTIM,  the  next  record  for  the  BGTFIL  file  (BGTTIM  and  BGTDEN)  is 
read  and  the  times  are  again  compared.  If  BGTTIM  exceeds  24.,  the  file  is 
rewound  before  it  is  read. 

REFERENCES: 

This  routine  is  called  from  COV,  the  main  routine,  early  in  the 
time  Iood. 
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3.3.3  Subroutine  BGTSCN 


PURPOSE: 

This  routine  Is  responsible  for  accumulating  the  quantity  of 
background  traffic  In  a  scan  area  bounded  by  points  A  and  B  along  a 
spacecraft  annulus  of  depth  DISTDP. 

INPUT  ARGUMENTS: 


The  following  Is  a  list  of  Input  arguments  passed  by  the  argument 

list: 


SCLAT 

SCION 

ALAT 

ALON 

BLAT 

BLON 

S 

L INNER 
01 STOP 


Spacecraft's  suborbital  point  latitude  In  degrees; 
Spacecraft's  suborbital  point  longitude  In  degrees; 
Latitude  of  one  scan  area  boundry  in  degrees; 

Longitude  of  one  scan  area  boundry  In  degrees; 

Latitude  of  other  scan  area  boundry  In  degrees; 

Longitude  of  other  scan  area  boundry  in  degrees; 

Earth  central  angle  of  annular  radius  in  degrees; 

Logical  variable,  TRUE  for  Inner  annulus,  FALSE  for  outer 
annul as  search; 

Depth  of  search  area  In  degrees 


OUTPUT  ARGUMENTS:  TOTTRF 


TOTTRF  Total  number  of  background  aircraft  traffic  found  in  the 


scan  area. 


METHOD: 


A  series  of  ten  degree  wide  subregions  Is  defined  within  the  scan 
area.  The  midpoints  of  these  subregions  and  their  areas  are  used  to 
define  the  density  blocks  for  routine  B6T.  The  final  subregion  may  be 
less  than  ten  degrees  In  width,  being  the  remainder  of  the  region.  The 
area  of  a  subregion  Is  defined  to  be  depth  times  the  width  times  the 
degree-to-nauti cal -mile  conversion  squared. 

REFERENCES: 

BGTSCN  Is  called  by  BRRIER  Subroutine. 

BGTSCN  calls  BGT,  BRNG,  DGLIM,  FLALOB  and  TRI. 


3.3.4  SUBROUTINE  BRRIER 


PURPOSE: 

Determines  If  a  satellite  detection  annulus  edge  intersects  an 
aircraft  barrier.  If  there  Is  an  intersection,  BRRIER  tests  to  see  if  any 
aircraft  are  detected  In  the  search  swath  along  the  annulus  edge  and  between 
the  two  points  of  Intersection.  Because  aircraft  barriers  can  be  annulur,  a 
satellite  detection  annulus  edge  can  Intersect  a  barrier  twice. 


Aircraft  barriers  are  searched  at  the  inner  and  outer  edge  of  the  annulus. 


INPUTS: 


BRREF 


Array:  See  BRBLD  for  description; 


TARGTS 

Array:  See  BRDLD  for  description; 

IBR 

Aircraft  barrier  ID; 

BELEVO 

Outer  annulus  edge  at  wheel  barriers  are  searched 
(degree  earth  grazing  angle); 

BELEVI 

Inner  edge  of  annulus  at  which  aircraft  barriers 
are  searched  (degrees,  earth  grazing  angle); 

TPOS  (1,NTAR) 

Aircraft  latitude; 

TPOS  (2,NTAR) 

Aircraft  longitude; 

TPOS  (4,NTAR) 

Time  since  aircraft  began  course; 

BW 

Beamwidth  of  radar  beam: 

OUTPUTS: 


NACS  Number  of  aircraft  barriers  being  searched  by 

current  spacecraft; 

NTSCH  Number  of  aircraft  and  ship  barriers  being 

searched  by  current  spacecraft;  * 

Locations  of  annulus  barrier  Intersections  defining  search  swaths  are 
stored  In: 

PT  (1,  NTSCH,  NSC)  Latitude  of  first  Intersection; 

PT  (2,  NTSCH,  NSC)  Longitude  of  first  Intersection; 

PT  (3,  NTSCH,  NSC)  Latitude  of  second  Intersection; 

PT  (4,  NTSCH,  NSC)  Longitude  of  second  Intersection; 

TPOS  (1,  BROW)  Latitude  of  first  Intersection; 

TPOS  (2,  BROW)  Longitude  of  first  Intersection; 

TPOS  (3,  BROW)  Latitude  of  second  Intersection; 

TPOS  (5,  BROW)  Longitude  of  second  Intersection; 

* 

Where  BROW  Is 

-  first  row  In  TARGTS  array  with  the  barrier  ID  (for  the  first  swath) 

-  second  row  In  TARGTS  array  with  this  barrier  ID  (for  the  second 
swath) 

Note:  TPOS  Is  changed  every  time  the  sensor  loop  Is  Incremented  while  PT  Is 
changed  only  when  the  time  loop  Is  Incremented.  PT  Is  used  to  write 
latitudes  and  longitudes  of  swath  limits  on  the  SCPLT  file. 

RES  (1,  NISCH)  Width  of  search  swath,  measured  In  spacecraft 

centered  degrees; 

RES  (2,  NISCH)  Search  cycle  time  (minutes); 

RES  (3,  NISCH)  Available  dwell  time  for  a  single  bean  (minutes); 

RES  (4,  NISCH)  Width  of  search  swath,  measured  In  rotational 

degrees  about  suborbital  point; 
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£2. 


RES  (5,  NISCH) 


RES  (6,  NISCH) 


Annulus  edge  bordering  search  swath  (earth 
central  angle  In  degrees); 

Number  of  aircraft  (Includes  background  traffic 
plus  threat  aircraft)  that  are  in  search  swath 
and  that  need  to  be  tracked. 


METHOD: 


BRRIER  performs  three  functions  for  each  aircraft  barrier  for  each 
spacecraft  sensor  at  each  time  step: 

1.  Determines  If  the  Inner  or  outer  edge  of  the  annulus  Intersects  the 
barrier;  thus  creating  a  search  swath; 

2.  Completes  the  entries  In  the  RES  array  for  each  search  swath.; 

3.  Determines  if  there  are  any  threat  aircraft  In  the  search  swath 
being  detected  for  the  first  time; 

If  there  are  two  search  swaths,  then  the  ID  for  the  first  one  Is  the 
odd  number  Input  as  the  barrier  ID  In  the  TAROTS  array  and  the  ID  of  the 
second  one  Is  the  Input  ID  number  plus  one. 

To  determine  If  there  Is  a  search  swath,  BRRIER  uses  CONES  to  test  for 
the  Intersections  of  the  annulus  edge  with  the  following  four  circles: 

1.  The  circle  defining  the  Inner  edge  of  the  barrier; 

2.  The  clrle  defining  the  outer  edge  of  the  barrier; 

3.  The  great  circle  defining  each  end  of  the  barrier.  See  description 
of  BRBLD. 

If  these  intersections  exist,  they  are  then  tested  against  the  barrier 
specifications  input  on  the  TARGTS  array.  This  Is  done  by  transforming  the 
Intersections  to  a  new  coordinate  system  whose  pole  Is  the  barrier  center  and 
whose  zero  meridian  Is  the  barrier's  bisector  (or  axis  of  symmetry).  After 
this  transformation,  the  colatitude  and  longitude  of  each  Intersection  are  the 
distance  from  the  barrier  center  and  the  bearing  relative  to  the  barrier 
bisector.  Then,  for  an  Intersection  to  be  on  the  barrier,  the  colatitude  must 
be  on  or  between  the  Inner  and  outer  edges  of  the  barrier  and  the  bearing  must 
be  less  than  or  equal  to  the  half-angle  of  the  barrier  width  (see  BRBLD). 
Intersections  passing  these  tests  are  Input  to  the  SRT  array  for  scheduling  by 
calling  FSTEP. 

The  entries  In  RES  are  then  computed  for  each  swath.  Search  swath 
width  In  rotational  degrees  Is  computed  by  converting  the  points  of 
intersection  to  a  new  coordinate  system  whose  pole  Is  the  spacecraft 
suborbital  point  and  whose  zero  meridian  Is  along  the  spacecraft's  Inertial 
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velocity  vector.  The  depth  of  the  search  swath  is  the  depth  of  the  radar  beam 
footprint  or  the  distance  traveled  by  the  suborbital  point  in  one  simulation 
time  step,  whichever  is  larger.  The  period  or  time  available  to  scan  (or 
search)  the  swath  is  the  time  it  takes  the  annulus  edge  to  travel  over  a 
distance  equal  to  the  depth  of  the  swath.  This  period  is  also  called  the 
search  cycle  time;  to  maintain  a  10W  probability  of  sweeping  over  an  aircraft 
as  it  goes  through  the  search  swath,  the  search  swath  must  be  scanned  at  least 
once  per  cycle  time.  Available  beam  dwell  time  is  the  ratio  of  cycle  time  to 
swath  width  (In  spacecraft  centered  degrees)  times  the  beamwidth.  The  number 
of  aircraft  tracked  Is  the  number  of  threat  aircraft  detected  (described 
below)  plus  one-half  the  number  of  background  aircraft  in  the  swath  (because 
It  Is  assumed  that  half  the  aircraft  in  a  barrier  are  headed  towards  the 

Soviet  side  of  the  barrier  and  half  are  headed  away  from  the  Soviet  side). 

Finally,  the  tests  to  determine  if  any  threat  aircraft,  which  are 
Input  In  the  TARGTS  array,  are  detected  in  the  swath  are: 

1.  Has  the  aircraft  begun  its  course?  If  yes,  then; 

2.  Has  the  aircraft  ended  its  course?  If  no,  then; 

3.  Is  the  aircraft  inside  the  search  swath?  If  yes,  then; 

4.  Does  the  aircraft  have  a  doppler  return  greater  than  the  minimum 
discernable  velocity? 

If  all  these  tests  are  passed,  then  the  aircraft's  position  relative 
to  the  spacecraft  is  input  to  the  SRT  array  by  calling  FSTEP. 
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•  3.3.5  SUBROUTINE  BRRIRS 

* 

PURPOSE: 

Models  detections  of  ship  barriers  by  radar  satellites.  A  ship  search 
swath  Is  defined  by  the  segment  of  the  outer  edge  of  an  annulus  that  Is  Inside 
a  ship  barrier.  Because  ship  barriers  are  always  circular,  there  Is  never 
more  than  one  search  swath  per  barrier. 

INPUTS: 

SHELEV  Annulus  edge  at  which  ship  searches  are  conducted. 

Expressed  In  terms  of  earth  grazing  angle  (degrees); 
IBR  Number  of  ship  barrier  being  tested; 

BREFS  ARRAY  See  BRBLD  for  description; 

NSC  Number  of  current  spacecraft; 

OUTPUTS: 

NSHS  Number  of  ship  barriers  being  searched  by  current 

spacecraft; 

NTSCH  Number  of  aircraft  and  ship  barriers  being  searched  by 

current  spacecraft; 

Each  ship  search  swath  Is  bounded  by  the  intersections  of  the  search  annulus 
and  the  circle  describing  the  barrier.  The  latitude  and  longitude  (In 
degrees)  of  each  Intersection  Is  stored  as  follows: 

TPOS  (1,  NTAR)  Latitude  of  first  Intersection; 

TPOS  (2,  NTAR)  Longitude  of  first  Intersection; 

TPOS  (3,  NTAR)  Latitude  of  second  Intersection; 

TPOS  (5,  NTAR)  Longitude  of  second  Intersection; 

PT  (1,  NTSCH,  NSC)  Latitude  of  first  Intersection; 

PT  (2,  NTSCH,  NSC)  Longitude  of  first  Intersection; 

PT  (3,  NTSCH,  NSC)  Latitude  of  second  intersection; 

PT  (4,  NTSCH,  NSC)  Longitude  of  second  Intersection; 


RES  (1,  NTSCH) 

RES  (2,  NTSCH) 
RES  (3,  NTSCH) 
RES  (4,  NTSCH) 

RES  (5,  NTSCH) 

RES  (6,  NTSCH) 


Width  of  search  swath,  measured  In  spacecraft 
centered  degrees; 

Search  cycle  time; 

Available  dwell  time  for  a  single  beam; 

Width  of  search  swath,  measured  In  rotational 
degrees  about  the  suborbital  point; 

Annulus  edge  bordering  search  swath;  expressed  as 
an  earth  central  angle  (In  degrees)  from  the 
suborbital  point; 

Not  used; 


ISCAN  (NTAR,  NSC) 
STBER  1,  STBER  1 
ST 


Set  to  one  If  barrier  is  searched  by  spacecraft; 
is  used  on  the  SCAN  output  file; 

Bearings  (in  degrees)  from  suborbital  point  to 
Intersections  of  annulus  edge  and  ship  barrier; 
Earth  central  angle  between  suborbital  point  and 
annulus  edge  at  which  barrier  search  Is 
performed; 


METHOD: 

Subroutine  CIRSCT  Is  called  to  compute  location  of  two  Intersections 
between  annulus  edge  at  which  ship  barriers  are  searched  and  circular  ship 
barrier.  If  there  Is  no  Intersection  then  CIRSCT  sets  IPOINT  to  zero,  which 
causes  BARRIRS  to  return  Immediately.  If  the  two  circles  do  Intersect,  then 

o  NSHS  and  NTSCHS  are  incremented ; 
o  The  latitude  and  longitude  of  the  Intersections  are 
stored  In  TPOS  and  PT; 

o  ISCAN  (NTAR,  NSC)  Is  set  to  one,  when  NTAR  Is  the 
row  of  the  TARGTS  array  describing  this  barrier 
and  NSC  Is  the  number  of  the  current  spacecraft; 
o  Subroutine  FSTEP  Is  called  to  enter  the  two  points 
of  Intersection  into  the  SRT  array  for  clustering 
and  scheduling; 

o  The  first  five  entries  In  RES  are  computed; 


3.3.6  SUBROUTINE  CLUSTR 
PURPOSE: 

Partition  all  the  search  areas  In  the  annulus  of  a  single  active 
satellite  Into  groups  (or  clusters)  which  can  be  mechanically  scanned  by 
moving  the  antenna  feed  of  the  satellite  boresight.  There  are  two  types 
of  search  areas,  barrier  search  areas  and  aircraft  reacquisition  and 
tracking  areas.  A  barrier  search  area  Is  a  swath,  one  beam  footprint 
deep,  along  the  edge  of  the  annulus  which  Is  Inside  a  barrier.  Aircraft 
reacquisition  and  tracking  areas  are  much  smaller  areas  of  uncertainty 
which  are  searched  to  detect  aircraft  that  have  been  previously 
detected.  For  each  cluster,  the  following  are  computed: 

o  boresight  cumpolnt 
o  background  traffic  In  search  areas 
o  energy  used  to  scan  cluster 


INPUTS: 

CLSWID  degrees  Angle,  centred  at  suborbital  point  and 

expressed  In  rotation  degrees  about  that 
point,  defining  the  maximum  width  of  a 
cluster,  l.e.,  the  maximum  width  that 
can  be  scanned  from  a  fixed  boresight 
position  by  moving  the  antenna  feed. 

NTD  Total  number  of  targets  detectble  in 

annulus.  In  the  context,  targets 
Include: 

o  Intersections  of  barrier 
boundaries  with  annulus 
o  aircraft  Inside  annulus 


AA  (35,35) 


NSC 

SRT  (3,35) 


A  square  symmetric  zero  diagonal  array 
containing  the  suborbital  point  centered 
angle  between  all  pairs  of  tarqets 
Inside  the  annulus.  For  barrier  search 
areas,  aircraft  not  yet  detected,  and 
ships  in  barrier  search  areas,  the 
diagonal  element  is  reset  from  zero  to 
the  ID  number  of  the  search  area. 

Number  of  current  spacecraft 

Stores  for  each  target  in  annulus: 

1.  Row  of  target  on  TARGTS  array 

2.  Bearing,  relative  to  spacecraft 
velocity  bearing,  from  suborbital 
point  to  target  (degrees) 

3.  Spacecraft  look  anqle,  from  nadir 
(degrees). 


OUTPUTS: 


CLM(36,10)  These  arrays  can  store  up  to  41  items 

CLM2(5,10)  of  Information  about  each  cluster: 

CLM3(25,10) 


CLM(1,ICL) 

CLK(2,ICL) 

CLM(3,ICL) 

CLM(4&5,ITCL) 


Cluster  number 

Number  of  spacecraft  viewing  cluster 
Cluster  number  for  this  spacecraft 
Bearings  from  spacecraft  suborbital 
point  to  targets  that  are  gratest  angle 
apart; 
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CLM(6&7,ITCL) 

CLM(8,ITCL) 

CLM(9,ITCL) 

CLM(10,ITCL) 

CLH(11,ITCL) 
CLM( 12-36, ITCL) 

CLM2(1, JCL) 
CLM2(2,JCL) 

CLM3(1,ITCL) 

CLM3(2,ITCL) 

CLM3(3,ITCL) 

CLM3(4,ITCL) 

CLM3(5,ITCL) 

CLM3(6,ITCL) 

CLM3(7,ITCl) 


Spacecraft  look  angles  to  see  targets 
that  are  greatest  angle  apart.  4-6  are 
for  the  target  tha  is  closest  to  the 
suborbital  point;  5-7  are  for  the  tarqet 
that  is  furthest  from  the  suborbital 
point; 

Bearing  to  boresight  AIM  point 
Spacecraft  look  angle  to  boresight  AIM 
point; 

Maximum  angular  separation  of  targets  in 
cluster; 

Number  of  targets  In  cluster 
Number  of  each  target  in  cluster; 
Initially  the  number  is  the  target's 
rank  in  the  SRT  array.  Maximum  number 
of  targets  currently  allowed  in  a 
cluster  is  25.; 

hours  Earliest  time  of  last  detection  across 

all  targets  in  the  cluster; 
hours  Minimum  duration  of  current  detection 

Interval  across  all  targets  In  cluster; 

Barrier  number; 

Earth  central  angle  of  search  (Inner  or 
outer  edge  of  annulus)  In  degrees; 

Search  cycle  time  (minutes)  from 
barrier; 

Search  amplitude  In  spacecraft  degrees; 
Search  time  (minutes),  hopefully  less 
than  cycle  time; 

Search  energy  (KJ)  expended  during  time 
step; 

Background  traffic  In  search  area  (Inc 
targets ) ; 

Time  P/l  on  during  search  (hrs); 


CLM3(8,ITCL) 


CLM3(9,ITCL) 

Number  of  targets  to  be  reacquired; 

CLM3(10,ITCL) 

Energy  (KJ)  used  to  reacquire  targets; 

CLM3(11,ITCL) 

Time  (minues)  to  reacquire  targets; 

CLM3(12,ITCL) 

Background  traffic  in  reacquisition 
areas  (inc  targets); 

CLM3(13,ITCL) 

Background  traffic  in  tracking  areas 
(inc  targets); 

CLM3(14,ITCL) 

Number  of  raid  targets  to  be  tracked; 

CLM3(15,ITCL) 

Track  time/cycle  time  ratio; 

CLM3(16,ITCL) 

Maximum  number  of  targets  that  can  be 

tracked  in  the  available  track  time; 

CLM3(17,ITCL) 

Track  energy  (KJ); 

CLM3(18,ITCL) 

Total  background  traffic  encountered; 

CLM3(19,ITCL) 

Time  P/L  on  during  tracking  (hrs); 

CLM3(20,ITCL) 

Time  P/L  on  during  reacquisition  (hrs); 

CLM3(21,ITCL) 

Total  energy  used  (KJ)  to  search, 
reacquire,  and  track  targets; 

CLM3(22,ITCL) 

Total  time  P/L  on  during  period; 

CLM3(23,ITCL) 

Not  used; 

CLM3(24,ITCL) 

Not  used; 

CLM3(25,ITCL) 

Not  used; 

CLCTR(2,5,8) 

CLCTR(1,IC2,NSC) 

degrees 

Longitude  of  boresight  aimoint; 

CLCTR(2,ICL,NSC) 

degrees 

Latitude  of  boresight  aimpoint; 

CLBG(5,5,9) 

degrees 

Stores  spacecraft  centered  angles 
between  all  pairs  of  boresight  aimpoints 
in  annulus; 

The  diagonal,  which  otherwise  would  be 
zero,  stores: 

CLBG(l(lf NSC) 

Is  #  If  there  is  a  scheduling  conflict, 
i.e.,  if  this  spacecraft  can  see  more 
than  one  cluster; 

CLBG(2,2,NSC) 

Number  of  clusters  in  annulus; 

CLBG(3,3,NSC)  Number  of  clusters  for  which  this 

spacecraft  Is  responsible.  (Another 
spacecraft  may  be  able,  for  example,  to 
see  one  of  the  clusters  visible  to  this 
spacecraft. ) 

Total  angle  to  visit  all  clusters. 


CLBG(4,4,N$C) 


METHOD: 


Subroutine  CLUSTR  partitions  th  targets  in  the  AA  array  into 
clusters,  which  are  defined  and  stored  in  the  CLM,  CLM2,  and  CLM3 
arrays.  Once  all  the  targets  In  a  cluster  are  identified,  the  boresight 
aimpoint  (in  the  center  of  the  cluster)  and  the  energy  required  to  scan 
the  cluster  are  computed. 

The  following  rules  guide  the  building  of  a  cluster: 


1.  A  cluster  contains  no  more  than  one  barrier  search  area  swath. 

2.  The  entire  swath  (defined  as  the  annulus  edge  between  two 
annulus  -  barrier  Intersections)  must  be  in  one  cluster,  even 
If  the  resulting  cluster  width  is  greater  than  CLSWID.  (If 
this  occurs,  the  spacecraft  must  slew  while  the  antenna  feed 
Is  moving). 

3.  A  cluster  will  generally  contain  an  aircraft  barrier  or  a  ship 
barrier  and  several  aircraft  track  or  reacquisition  search 
areas  which  are  approximated  by  the  current,  true  position  of 
the  aircraft. 

4.  Clusters  are  first  built  around  barrier  search  swaths  and  then 
around  unclustered  aircraft.  The  first  member  of  a  cluster  is 
called  the  pivot. 

5.  If  the  pivot  Is  a  barrier  search  swath  Intersection  or  as  yet 
undetected  aircraft  or  ship  appearing  in  a  barrier  search 
swath,  then  all  other  members  of  that  barrier  search  swath  are 
included  In  the  cluster.  CLSWID  is  replaced  by  the  width  of 
the  barrier  search  swath  if  it  Is  greater  than  the  search 
swath. 

6.  After  the  pivot  target(s)  are  in  the  cluster,  all  other 
aircraft  in  the  AA  array  are  tested  to  see  If  they  can  also  be 
in  the  cluster.  To  be  in  the  cluster,  the  candidate  aircraft 
must  be  less  than  CLSWID  degrees  away  from  all  targets  already 
In  the  cluster. 

7.  The  process  of  identifying  a  pivot  and  testing  all  other 
unassigned  AA  array  entries  against  that  pivot  continues  until 
all  entries  of  the  AA  array  are  clustered. 

8.  The  diagonal  element  In  the  AA  array  for  an  individual  target 
is  set  to  99  when  that  target  has  been  Included  in  a  cluster. 


|  The  boreslght  Is  aimed  at  a  spacecraft  look  angle  that  bisects  the 

spacecraft  centered  angle  from  the  Inner  to  the  outer  edge  of  the  annulus 
and  at  bearing  angle  that  bisects  the  width  of  the  cluster. 

The  logic  to  compute  background  traffic  In  search  areas  and  the 
time  and  energy  to  scan  a  cluster  Is  still  developmental  and  Is  not 
correct  as  of  this  writing. 

* 

i 

\ 


I 
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I 
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3.3.7  SUBROUTINE  CLSTPP 
PURPOSE:  * 

Subroutine  CLSTPP  (cluster  port  processor)  Is  the  scheduler  for  active 
satellites.  If  an  active  satellite  can  view  more  than  one  cluster  (see 
description  of  subroutine  CLUSTR),  then  the  schedule  decides  at  which  cluster 
to  aim  the  boreslght. 

INPUTS: 

CLBG  (5,5,9)  See  description  of  subroutine  CLUSTR.  CLSTPP 

checks  to  see  If  two  satellites  see  the  same 
cluster.  In  which  case  the  third  diagonal  element 
Is  decremented  for  one  of  the  spacecraft  and  the 
first  diagonal  element  may  be  set  to  a  blantk. 

CLM  (36,10)  See  description  in  CLUSTR.  For  each  cluster 

(l.e.,  for  each  row  on  the  CLM  array),  CLSTPP 
uses: 

CLM  (1,ICL):  Cluster  number; 

CLM  (2,ICL):  Number  of  spacecraft 

viewing  cluster; 

CLM  (3,ICL):  Cluster  number  for  this 
spacecraft; 

CLM  (11,ICL):  Number  of  targets  in 
this  cluster; 

CLM  (12-36):  Target  number  (l.e.,  row 
on  targets  array)  of 
each  target  In  cluster; 

Note:  In  this  context,  "target"  can  mean  the  endpoints  of  a  barrier  search 
swath. 


CLM3  ( 1 , ICL) 


Earliest  time  of  last  detection  across  all 
targets  In  cluster. 


CLM2  (2,ICL) 

CLM3  (1,ICL) 
XLN  (NTAR.l) 
XLN  (NTAR.ll) 


NSC  I 

OUTPUTS: 

SCSCH  (9,8) 

SCSCH  (l.NSL) 

SCSCH  (2, NSC) 

SCSCH  (3, NSC) 

SCSCH  (4, NSC) 

SCSCH  (5, NSC) 

SCSCH  (6  and  7.NSC) 

SCSCH  (8, NSC) 

SCSCH  (9, NSC) 

ICLST  (9,  10) 


Minimum  duration  of  current  detection  Interval 
across  all  targets  In  cluster. 

10  number  of  barrier  search  swath 
Time  (hours)  at  which  last  detection  occurred. 
Time  (hours)  at  which  current  detection  began;  Is 
set  to  zero  If  target  was  not  detected  at  last 
time  step. 

Total  number  of  clusters  seen  by  all  active 
sensors  (e.g.  active  spacecraft,  0T4  B  radars, 
and  DEW  line)  at  current  time  step.  (By 
definition,  there  Is  no  more  than  one  cluster  per 
fixed  sensor.) 

The  number  of  active  sensors. 


Shows  how  each  active  sensor  Is  scheduled. 

1  If  sensor  sees  no  more  than  one  cluster; 

0  If  sensor  sees  none  or  several  clusters; 

1  If  sensor  has  been  assigned  to  a  cluster; 

0  If  sensor  has  not  been  assigned  to  a  cluster; 
Row  In  CLM/CLM2/CLM3  arrays  containing  cluster  to 
which  this  sensor  has  been  assigned; 

Number  of  cluster  to  which  spacecraft  has  been 
assigned.  This  number  can  be  less  than  SCSCH 
(3, NSC),  but  only  If  two  or  more  sensors  see  the 
same  cluster  (l.e.,  have  matching  clusters). 
Number  of  targets  In  cluster; 

First  and  last  rows  In  CLM/CLM2/CLM3  containing 
clusters  that  can  be  seen  by  this  spacecraft; 
Earliest. time  of  last  detection  across  all 
targets  In  assigned  cluster; 

Minimum  duration  of  current  detection  Interval 
across  all  targets  In  assigned  cluster; 

Shows  which  active  sensors  are  assigned  to  each 
cluster; 
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ICLST  (1,ICL)  Total  number  of  active  sensors  assigned  to  this 
cluster; 

ICLST  (2-9, ICL)  Active  sensor  ID  numbers  of  sensors  scheduled  to 
view  this  cluster; 

1CLW1D  (8)  Width  of  cluster  (In  terms  of  rotation  degrees 

about  spacecraft  suborbital  point)  which  has  been 
scheduled  for  each  spacecraft; 

METHOO: 

CLSTPP  assigns  or  schedules  active  sensors  to  look  at  one  cluster  per 
time  step.  For  OTH  B  radar  and  DEW  Line  this  Is  trivial,  because  by 
definition  OTH  B  and  the  DEW  Line  each  see  only  one  cluster.  Scheduling 
spacecraft  Is  harder,  because.  In  general,  each  spacecraft  has  at  each  time 
step  more  than  one  cluster  (In  Its  annulus)  that  needs  to  be  scanned.  Recall 
that  a  cluster  Is  a  group  of  aircraft  and  up  to  one  aircraft  or  ship  barrier 
search  swath  that  can  be  scanned  from  a  fixed  boreslght  position  by  moving  the 
antenna  feed.  When  an  active  spacecraft  has  a  scheduling  conflict,  the 
conflict  Is  resolved  based  on  the  following  priorities  (listed  In  descending 
order  of  Importance). 

1.  Do  not  sche&le  a  cluster  If  It  Is  already  scheduled  to  another 
sensor  (Implies  that  two  sensor  see  the  same  cluster)  unless  the 
current  sensor  has  In  Its  cluster  a  target  which: 

-  Is  not  seen  by  the  other  sensor; 

•  has  longer  detection  gap  than  the  targets  by  the  other  sensor; 

-  has  a  shorter  current  detection  period  than  that  of  any  of  the 
targets  seen  by  the  other  cluster; 

2.  Choose  the  cluster  that  has  target  with  the  longests  detection  gap 
unless  the  maximum  gap  Is  less  than  .9  hrs. 

3.  Choose  the  cluster  which  has  the  target  with  the  shortest  detection 
period  If: 

-  the  target  was  observed  at  the  last  time  step; 

>  The  detection  period  Is  less  than  .3  hrs; 


4.  Choose  the  cluster  that  has  the  most  targets. 

It  occasionally  happens  that  a  sensor  Is  not  scheduled  to  see  a 
cluster  -  all  the  clusters  It  sees  are  already  scheduled  to  other  sensors. 

The  scheduler  assigns  sensors  In  ascending  order  -1,  2,  ...»  NSCI;  It  does  not 
unassign  nor  reassign  sensors.  (This  would  be  a  desirable  Improvement.)  When 
a  sensor  Is  scheduled  to  a  particular  cluster,  this  Is  recorded  In  the  SCSCH 
array. 


3.3.8  Subroutine  ECLPSF 
PURPOSE: 

This  routine  Is  responsible  for  updating  the  SUNVEC  data  to  the 
current  simulation  date/tlme.  It  Is  called  within  the  early  stages  of  the 
Time  Loop  of  COV.  The  SUNVEC  Is  used  to  determine  the  eclipse  status  of 
the  spacecraft  within  that  Time  Loop. 

INPUT  ARGUMENTS: 

DATE  The  current  Julian  date  and  time,  the  latter  expressed  as 
a  fraction  of  a  day. 

OUTPUT  ARGUMENTS: 

SUNVEC  The  x-,  y-  and  z-  unit  vector  for  the  sun's  position 
relative  to  the  Earth. 

METHOD: 

The  algorithm  Is  embeded  In  SUN  and  the  other  utility  routines 

called. 

REFERENCES: 

ECLPSF  Is  called  by  COV.  It  calls  the  following  utility 
routines:  SUN,  OTOC,  EQEC,  and  RMMUL. 


’ZfrZxk  %. 
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3.3.9  Function  ECLPSF 
PURPOSE: 


ECLPSF  Is  a  logical  function  Intended  to  return  a  TRUE  value  when 
the  spacecraft  position  Input  Is  In  solar  eclipse  during  the  current  time, 
and  otherwise,  FALSE. 

INPUT  ARGUMENTS: 


SUNVEC  The  x-,  y-  and  z-  unit  vector  of  the  sun's  position 
relative  to  the  Earth  at  the  current  time. 

SCLON  Longitude  of  suborbital  point  in  degrees 
SCLAT  Latitude  of  suborbital  point  In  degrees 
SCRAD  Orbital  radius  In  feet. 

OUTPUT  ARGUMENTS: 

ECLIPSE  A  logical  variable,  set  TRUE  for  eclipse,  otherwise,  V 

FALSE. 


METHOD: 

The  unit  vector  for  the  spacecraft  Is  computed  and  the  dot  product 
of  this  vector  with  the  solar  unit  vector  Is  taken.  From  the  comparison 
of  the  result  with  the  horizon  of  the  spacecraft,  the  eclipse  state  can  be 
determined. 

REFERENCES: 

The  routine  Is  used  in  COV.  It  calls  DOT  and  MAGN. 
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3.3.10  SUBROUTINE  FSTEP 


PURPOSE: 

Performs  the  first  step  of  active  satellite  clustering  and  scheduling, 
which  Is  to  store  the  bearing  and  elevation  angles  of  detectable  search  areas 
in  the  SRT  array  in  order  of  their  bearing  relative  to  the  spacecraft's 
velocity  heading. 

INPUTS:  (all  angles  In  degrees) 

Earth  central  angle  from  suborbital  point  to  a  search 
area  location; 

Bearing  (from  north)  from  suborbital  point  to  search 
area  location; 

Spacecraft  velocity  bearing; 

Search  area  10  number  (row  on  TARGTS  array  that 
describes  the  search  area).  A  search  area  can  be  a 
swath  through  an  aircraft  or  ship  barrier  or  the 
estimated  current  position  of  an  aircraft. 

OUTPUTS:  (all  angles  In  degrees) 

SRT ( 1 ,  I)  Search  area  ID  number; 

SRT ( 2 ,  I)  Bearing,  relative  to  spacecraft  velocity  heading,  from 

spacecraft  suborbital  point  to  search  area  location; 
SRT(3,  I)  Spacecraft  look  angle  from  nadir  to  target; 

Note:  all  detectable  search  areas  are  entered  In  the  SRT  array  In  ascending 
order  based  on  bearing  from  suborbital  point  to  search  area.  Bearings  can 
vary  from  -180°  to  +  180°. 

NTO  Number  of  search  areas  detectable  by  current 

spacecraft; 


ST 

STBER 

SVBER 

NTAR 


METHOO: 


This  subroutine  Is  an  Intermediate  step  which  Is  prerequisite  for 
subroutine  CLUSTR.  Subroutine  CLUSTR  groups  detectable  search  areas  for  an 
active  satellite  into  groups  (called  clusters)  which  can  be  scanned  from  a 
fixed  boresight  position  by  moving  the  antenna  feed.  To  do  this,  CLUSTR 
requires  that  all  detectable  search  areas  be  stored  In  clockwise  order  around 
the  annulus.  FSTEP  puts  search  areas  in  this  order.  A  search  area  is: 

o  An  aircraft  barrier  search  swath; 

o  A  ship  barrier  search  swath; 

o  The  current  position  of  an  aircraft  that  has  already 
been  initially  detected; 

For  each  search  swath,  two  points  are  stored  In  SRT,  namely,  the  two 
Intersections  of  the  barrier  with  the  annulus  edge  bordering  the  search 
swath.  Currently,  the  model  approximates  a  search  area  to  reacquire  a 
previously  detected  aircraft  by  the  aircraft's  current,  actual  position.  This 
approximation  is  not  valid  If  the  aircraft  maneuvers  between  detections;  the 
search  area  should  be  based  on  an  a  position  extrapolated  from  the  last 
observed  position  and  should  have  some  area  of  uncertainty. 

A  search  area  Is  detectable  If  it  Is  In  the  spacecraft's  detection 


annulus. 


3.3.11  SUBROUTINE  FXSENS 


PURPOSE:  | 

Models  detections  of  aircraft  by  fixed  sensors,  e.g.,  Over-the-Horizon 
Backscatter  (OTH  B)  radar  and  the  DEW  Line.  The  approach  is  to  model  the 
areas  covered  by  these  two  systems  as  annular  sectors  and  to  compute  when 
aircraft  are  in  these  sectors. 

INPUTS: 

(All  angles  in  degrees;) 


SLAT 

SLON 

SHO 

SWD 

SIR 

SOR 


Latitude  of  annulus  center; 

Longitude  of  annulus  center; 

Bearing  from  north  of  axis  of  symmetry  (or 
bisector)  of  annular  sector; 

Angulur  width  of  annular  sector  (half  angle); 
Distance  from  annulus  center  to  Inner  edge  of 
annular  sector  (n.ml.); 

Distance  from  annulus  center  to  outer  edge  of 
annulus  sector; 


TARGTS  (NTAR,  1)  Greater  than  zero  if  target  is  an  aircraft; 

TARGTS  (NTAR,  7)  Time  at  which  aircraft  begins  flight; 


TP0S(1,  NTAR) 
TP0S(2,  NTAR) 
TPOS(3,  NTAR) 
TP0S(4,  NTAR) 


Aircraft's  current  latitude; 
Aircraft's  current  longitude; 
Aircraft’s  current  bearing; 

Time  since  aircraft  started  course; 


OUTPUTS: 

CLM( 1 ,  J) 
CLM(2,  J) 
CLM(3,  J) 

CLM( 1 1 ,  J) 


Cluster  number  (usually  equal  to  J) ; 

Sensor  number; 

Number  of  clusters  seen  by  this  sensor  (Is  always 
one  for  OTH  B  radar  and  DEW  Line); 

Number  of  aircraft  seen  by  this  sensor; 
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CLM  (12-35,  J)  Aircraft  ID  (l.e.,  row  In  TARGTS  array)  for  each 

aircraft  seen  by  sensor; 

ME^Hof: 

To  determine  If  an  aircraft  Is  In  the  coverage  area  of  an  OTH  B  radara 
or  the  DEW  Line,  FXSENS  checks  all  rows  In  the  TARGTS  array  to  find  if: 

o  row  represents  an  aircraft;  if  yes,  then; 
o  has  the  aircraft  begun  its  flight?  if  yes,  then; 
o  has  the  aircraft  ended  Its  flight;  If  not,  then  the 
current  position  of  the  aircraft  Is  transformed 
to  a  new  coordinate  system  whose  pole  Is  the 
center  of  the  annular  coverage  sector  (e.g.  the 
OTH  B  radar  transmitter  Is  the  center  of  the 
annular  sector  covered  by  the  OTH  B  radar  system) 
and  whose  zero  meridian  Is  the  axis  of  symnetry 
of  the  annular  sector.  In  the  new  coordinate 
system,  the  aircraft's  colatitude  Is  Its  distance 
■%  from  the  center  and  Its  longitude  Is  Its  angulur 

separation  from  the  center  of  the  annular 
coverage  sector.  The  tests  to  determine  If  the 
aircraft  Is  In  the  annular  coverage  sector  are, 
therefore: 

o  Is  aircraft's  new  colatitude  greater  than  the 
distance  to  the  Inner  edge  of  the  coverage 
sector;  If  yes,  then; 

o  Is  the  new  colatitude  less  than  the  distance  to  the 
outer  edge  of  the  coverage  sector;  If  yes,  then; 
o  Is  the  absolute  value  of  the  aircraft's  new 

longitude  less  than  the  ha if -angle  width  of  the 
coverage  sector; 

If  all  these  tests  are  passed,  then  the  aircraft  Is  In  the  coverage 
area  of  the  current  fixed  sensor. 


If  the  fixed  sensor  sees  one  or  more  aircraft,  then  FXSENS  creates  a 
cluster,  similar  to  the  clusters  created  for  spacecraft  by  subroutine 
CLUSTR..  However,  whereas  for  spacecraft  there  may  be  several  clusters  of 
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which  only  one  can  be  observed  at  a  time;  for  OTH  B  radar  or  the  DEW  Line, 
there  is  only  one  cluster,  i.e.,  all  aircraft  in  the  coverage  area  can  be 
seen  in  one  simulator  time  step.  A  cluster  is  created  by  making  the  entries 
shown  under  outputs  in  a  new  row  of  the  CLM  array.  This  makes  it  unnecessary 
to  call  CLUSTR  for  fixed  sensors.  However,  subroutine  CLSTPP,  which  schedules 
spacecraft  to  look  at  the  highest  priority  cluster,  uses  the  CLM  array  to 
assign  spacecraft  to  clusters  which  do  not  have  aircraft  that  are  seen  by  a 
fixed  sensor.  Thus,  fixed  sensors  can  simplify  spacecraft  scheduling. 


3.3.12  SUBROUTINE  JMEQ 


PURPOSE: 

Computes  an  aircraft's  velocity  component  along  the  line  of  sight  from 
the  aircraft  tothe  spacecraft. 

INPUTS:  (all  angles  in  degrees) 

ST  Earth  central  angle  from  suborbital  point  to  aircraft 

position; 

DBER  Suborbital-point-centered  angle  between  spacecraft 

velocity  bearing  and  great  circle  to  aircraft 
position; 

VT  Aircraft  velocity  (knots) 

OUTPUTS: 

DDOP  Aircraft's  velocity  component  along  the  line  of  sight 

from  the  aircraft  to  the  spacecraft  (knots); 

METHOD: 

Compute  the  aircraft  centered  angle  between  the  aircraft  velocity 
vector  and  the  aircraft-to-spacecraft  line-of-sight  radial.  This  is  done 
using  a  right  spherical  triangle  on  a  unit  sphere  whose  center  is  the 
aircraft.  Two  adjacent  (and  perpendicular)  sides  are  the  arcs  corresponding 
to  DBER  and  the  elevation  angle  (above  the  aircraft's  local  horizon)  to  the 
spacecraft.  The  third  arc  length  corresponds  to  the  angle  between  the 
velocity  vector  and  the  line-of-siqht  from  the  aircraft  to  the  spacecraft. 

The  aircraft  velocity  component  along  the  line-of-sight  is  its  speed  times  the 
cosine  of  the  angle  just  computed. 
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3.3.13 


SUBROUTINE  PWRBAL 


PURPOSE: 

This  routine  Is  responsible  for  maintaining  the  state  of  charge  of 
the  battery.  The  state  of  charge  is  updated  each  time  increment  based 
upon  the  amount  of  time  the  payload  Is  used  during  the  period  and  the 
solar  eclipse  state. 

INPUT  ARGUMENTS: 

The  following  parameters  are  input  via  the  calling  sequence: 

TUSED  Amount  of  time  payload  was  used  in  hours. 

ECLIPS  Logical  variable  indicating  eclipse  state,  TRUE  for 
eclipse;  FALSE,  otherwise. 

SOC  State  of  charge,  as  a  fraction,  at  start  of  time 
Interval. 

The  following  parameters  are  input  via  Common  X2  or  X8 

BATCAP  Maximum  battery  capacity  in  kilowatt  hours. 
r-  SOLAR  Power  input  from  solar  panels  in  kilowatts 
PPLON  Power  consumed  while  payload  is  on,  including  power 
consumption  of  platform,  in  kilowatts. 

PPLOFF  Power  consumed  while  payload  Is  off.  Including  power 
consumption  of  platform,  in  kilowatts. 

DUTCYC  Duty  cycle  of  payload. 

OTIME  Time  interval  In  hours. 

OUTPUT  ARGUMENTS: 

SOC  The  state  of  charge  at  the  end  of  the  time  interval 


METHOD: 


The  input  power  Is  determined  by  the  eclipse  state,  0  for  TRUE, 
PSOLAR  for  FALSE.  If  the  payload  was  used  during  the  period,  the 
consumption  rate  is  set  to  PPLON  for  a  period  of  DTYCYC  *  TUSED.  The 
period  the  payload  was  off  is  either  the  entire  interval  or  the  remainder 
of  the  interval  (DTIME-DTYCYC  *  TUSED)  and  power  consumption,  PPLOFF.  The 
net  power  is  then  computed  (Power  in  minus  Power  out)  as  well  as  an  ETA 
factor.  If  the  battery  is  discharging,  the  change  in  state  of  charge  is 
the  net  change  times  the  period  divided  by  ETA  times  battery  capacity.  If 
the  battery  Is  charging,  the  change  In  state  of  charge  is  equal  to  net 
change  times  ETA  times  .7576  times  the  period  divided  by  the  battery 
capacity.  The  new  SOC  Is  then  the  old  SOC  plus  the  change  in  state  of 
charge.  If  the  payload  was  used,  there  Is  a  period  that  it  was  not  used 
during  the  interval.  The  power  consumption  during  that  period  is  computed 
next  as  already  described  for  the  payload  off. 

REFERENCES: 


This  routine  is  called  by  COV 


3.4  OUTPUT 


3.4.1  SUBROUTINE  BGNCRS 
PURPOSE: 

BGNCRS  calls  SETL  EPSl  before  moving  target  begins  course  so  that 
coverage  statistics  computed  by  subroutine  summary  and  the  post-processor 
STATCV  Include  gaps  and  detections  occurring  at  the  beginning  of  the 
target's  course.  (Note  that.  In  the  statistics,  the  Initial  gap  or 
detection  will  be  EPSl  plus  the  duration  of  the  gap  or  detection.) 

ENDCRS  calls  SETL  EPSl  after  moving  target  begins  course  so  that 
coverage  statistics  computed  by  subroutine  summary  and  the  post-processor 
STATCV  Include  gaps  and  detections  occurring  at  the  end  of  the  target's 
course.  (Note  that.  In  the  statistics, the  final  gap  or  detection  will  be 
EPSl  dIus  the  duration  of  the  gap  or  detection.) 

INPUTS: 


STPROF 

hours 

Time  (since  start  of  simulation)  that  moving 
target  begins  course. 

EDPROF 

hours 

Duration  of  target's  course. 

EPSl 

hours 

See  above.  Normally  stt  to  .0001  hours. 

NSMCON 

- 

Number  of  target -ground  station  sets  on 
ISMTAB  array. 

NTRGTT 

- 

Number  of  targets  on  the  TARGTS/ITRGTS  array 
plus  NSMCON. 

DTIME 

hours 

Time  steo  size  for  simulation. 

TIME 

hours 

Current  time  being  simulated. 

OUTPUTS: 

Added  to  TIME;  Is  the  time  at  which  SETL  Is 
called  and  results  In  SETL  being  called  EPSl 
before  course  begins  or  EPSl  after  course 
ends.  Is  reset  to  999  at  end  of  subroutine. 


TAOJ 


hours 


METHOD: 


The  main  program,  COV,  calls  BGNCRS/ENDCRS  two  time  steps 
before/after  course  begins/ends.  For  BGNCRS,  a  modular  divide  determines 
how  much  time  elapses  In  the  next  time  step  before  the  course  begins;  EPS1 
Is  subtracted  from  the  time  period  to  compute  TADJ.  For  ENDCRS,  a  modular 
divide  determines  how  much  time  elapsed  In  the  previous  time  step  before 
the  course  ended;  EPS1  Is  added  to  this  period  to  compute  TADJ.  After 
SETL  Is  called,  TADJ  Is  reset  to  999. 

REFERENCES: 


Calling  routines:  COV 
Called  routines:  SETL 


3.4.2  Subroutine  CTOALF 


PURPOSE: 

This  routine  is  responsible  for  converting  the  numeric  cluster 
number  found  in  ISCAN  array  to  either  an  alpha  or  a  digit  (in  alpha 
form).  The  SCSCH(4,  I)  volumn  Is  also  converted  to  alpha.  The  ISA  vector 
is  set  from  the  rounded  value  of  the  CLBG (5, 5, I)  vector.  These  values  are 
used  in  the  SCAN  printout. 

INPUT  ARGUMENTS: 

The  following  Inputs  are  passed  via  the  calling  sequence: 

NSC1  Number  of  aopropriate  S/C; 

NTRGTS  Number  of  targets; 

SCSCH  Array  whose  (4,  I)  column  represents  the  selected 
cluster  id; 

The  ISCAN  array  and  the  CLBG  array  are  passed  via  the  X4  Common. 
OUTPUT  ARGUMENTS: 

The  following  outputs  are  passed  via  the  calling  sequence: 

I SCSCH  Array  whose  (4,  I)  column  represents  the  converted  ID  for 

the  selected  cluster; 

ISA  Vector  containing  the  rounded  value  of  the  CLBG  (5,  5,  I) 
column. 

The  ISCAN  array  is  also  used  for  output. 


METHOD: 


ISA  Is  derived  from  the  fol lowing: 

ISA(I)  -  IFIX(.5  +  CLBG(5,  5,  I)) 

ISCAN  Is  derived  from  the  following: 

If  ISCAN  (J,I)  -  0 

-ISCAN(J,I)  as  an  alpha  coded  digit  (1-9) 
ISCAN(J.I)  If  ISCAN(J,I)  <  0 

•A'-'Z*  If  0  <  ISCAN(J.I)  <  26 
If  ISCAN(J.I)  >  26 

ISCSCH(4,I)  Is  transformed  In  the  same  way  ISCAN  Is. 
REFERENCES: 

This  routine  Is  called  from  COV. 


3.4.3  SUBROUTINE  SETL 


PURPOSE: 

Gathers  statistics  on  the  duration  of  detections  and  of  gaps  between 
detections.  At  the  start  of  each  new  detection,  SETL  also  writes  on  unit  8 
the  time  at  which  the  detection  begins  and  the  duration  of  the  coverage  gap 
since  the  last  detection. 

INPUTS: 

\ 

\ 

TfNE  Current  clock  time  (hours); 

DTIMfc  Simulator  time  step  (hours); 

DTIMX  Larger  time  step  if  two  types  of  spacecraft  are 

modeled; 

SATTAK  Time  at  which  current  replication  started  (hours); 

TADJ  Amount  by  which  time  has  been  adjusted  by  ENCRS  or 

BGNCRS; 

NTAR  I  (row  on  TARGTS  array)  for  object  being  detected; 

IPL0T2  Must  be  equal  to  or  greater  than  one  to  cause  a  write 

on  unit  8; 

OELTIM  Time  since  aircraft  began  flight  or  since  start  of 

replication  (for  fixed  targets); 

NUMCT  Number  of  simultaneous  contacts  required  for 

detection; 

OUTPUTS:  (all  times  In  hours) 

XLN(I,  1)  Time  at  which  last  contact  occurred; 

XLN(I,  2)  Maximum  time  out  of  contact; 

XLN(I,  3)  End  time  of  maximum  time  out  of  contact; 

XLN(I,  4)  Maximum  time  out  of  contact; 

XLN(I,  5)  End  of  minimum  time  out  of  contact; 

XLN(I,  6)  Total  time  out  of  contact; 

XLN(I,  7)  Number  of  contacts; 

XLN(I,  8)  Sum  of  squares  of  down  times; 

XLN(I,  9)  Time  at  which  first  detection  began; 


XLN( I ,  10)  Time  at  which  last  detection  began; 

XLN(I,  11)  Tine  at  which  current  contact  began  (zero  if  currently 

out  of  contact); 

Data  written  on  unit  8: 

NTAR 
IATTAK 
DLTIM 
DT 

ICLS  (NTAR) 

METHOD: 

SETL  Is  called  everytlme  a  ground  station,  an  aircraft,  or  a  barrier 
Is  detected.  If  the  detection  Is  the  first  detection  after  a  gap,  then  the 
SLM  array  Is  updated  and  unit  8  written  on;  otherwise  a  RETURN  Is  executed. 
Special  tests  are  required  If  a  detection  occurs  exactly  at  the  start  of  a 
replication  or  at  the  first  time  step  Into  a  replication.  BGNCRS  and  ENDCRS 
call  SETL  just  before  and  after  aircraft  flights  so  gaps  before  the  first 
detection  and  after  the  last  detection  are  counted  In  the  statistics. 

However,  when  ENDCRS  calls  SETL,  DT  Is  always  less  than  the  time  step,  which 
would  cause  an  Immediate  return.  However,  this  return  Is  prevented  If  TADJ  Is 
negative,  and  TADJ  Is  negative  only  If  SETL  Is  called  from  ENDCRS. 

Special  format  statements  are  used  to  write  a  zero  before  single  digit 
positive  numbers  so  sort  utilities  will  sort  on  a  zero  Instead  of  a  blank. 


ID  of  aircraft  or  barrier  being  detected; 
Replication  number; 

DELTIM  (see  above); 

Time  since  last  contact  ended; 

1  If  aircraft  Is  classified;  0  If  otherwise; 


3.4.4  Subroutine  STATPT 


PURPOSE: 

This  routine  calls  STHIST  for  each  of  five  terms  to  printout  the 
post-run  statistics.  The  terms  are  abbreviated  SA,  CW,  SC,  EB,  and  RC. 

INPUT  ARGUMENTS: 

The  STATS  Common  provides  all  the  inputs  to  this  routine.  For  each 
of  the  five  terms,  the  following  data  are  provided  where  **  is  an 
abbreviation  for  one  of  the  five  terms: 


NCNT** 

Total  of  all  frequency  counts 

**MIN 

Minimum  value 

**MAX 

Maximum  value 

**SUM 

Sum  of  values 

**SUM2 

Sum  of  values  squared 

I**VEC 

Frequency  count  vector 

# 

OUTPUT  ARGUMENTS: 

For  each  term,  the  following  data  are  returned: 

**SUM  Mean 

**SUM2  Stand  deviation 


METHOD: 

The  routine  contains  data  statements  for  the  labels  for  each  of  the 
five  terms.  The  routine  calls  STHIST  for  each  term  after  heading  each 
term's  report. 

REFERENCES: 

STATPT  Is  called  by  COV.  It  calls  STHIST. 
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3.4.5  Subroutine  STATUD 


PURPOSE: 

This  routine  Is  responsible  for  maintaining  the  statistics  for  each 
of  the  five  terms.  The  count,  minimum,  maximum,  sum  and  sum  square  are 
maintained  and  the  frequency  bins  are  Incremented  as  appropriate.  The 
terms  are  abbreviated  SA,  CW,  SC,  EB  and  RC. 

INPUT  ARGUMENTS: 

The  calling  sequence  provides  the  fallowing  arguments: 

MODE  The  term  Identifier  (1  for  SA,  2  for  CW,  3  for  SC,  4  for 
EB  and  5  for  RC) 

VALUE  The  value  of  the  term. 

The  STATS  Common  Is  used  to  hold  these  statistics,  as  defined  in 

STATUP. 

OUTPUT  ARGUMENTS: 

The  updated  STATS  Common  Is  the  only  output  for  this  routine. 


METHOD: 

Depending  upon  the  value  of  MODE,  the  statistics  are  updated  for 
the  term.  For  the  basic  statistics,  the  formulae  are  as  follows  where  ** 
represents  the  term's  abbreviation: 

NCNT**  *  NCNT**  +  1 
**MIN  =  AMIN1(**MIN,  VALUE) 

**MAX  *  AMAX1(**MAX,  VALUE) 

**SUM  «  **SUM  +  VALUE 
**SUM2  »  **SUM2  +  VALUE  *  VALUE 
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For  this  frequency  array,  the  maintenance  is  dependent  upon  MODE. 
For  MODE  =  1, 

I  =  MIN0(18,  MAX0(1,  1  +  IFIX ( -1 . E-10  +  VALUE/10.))) 

For  MODE  =  2, 

I  *  MIN0(11,  MAX0(1,  1  +  I  FIX (-1 . E-10  +  VALUE/10.))) 

For  MODE  =  3, 

I  =  MINO(12,  MAX0(2,  2  +  IFIX(-1.E-10  +  VALUE/10.))) 
or  I  =  1  if  VALUE  <  0 
For  MODE  *  4, 

I  =  MIN0(11,  MAX0{1,  1  +  IFIX (-1 . E-10  +  VALUE/10.))) 

For  MODE  =  5, 

I  =  MIN0(9,  MAXO(l,  1  +  IFIX(-1.E-10  +  VALUE/25.))) 
and  I***VEC  (I)  *  I***VEC(I)  +  1 

REFERENCES: 


The  routine  is  called  by  COV. 


3.4.6  SUBROUTINE  SMMRY 


PURPOSE: 


Computes  statistics  summarizing  one  replication  using  data  stored  in 
the  XLN  array. 

INPUTS 


OUTPUTS 


XLN ( I „  1)  Time  at  which  last  contact  occurred; 

XLN( 1 ,  2)  Maximum  time  out  of.  contact; 

XLN(I,  3)  End  time  of  maximum  time  out  of  contact; 

XLN(If  4)  Maximum  time  out  of  contact; 

XLN(I,  5)  End  of  minimum  time  out  of  contact; 

XLN(I,  6)  Total  time  out  of  contact; 

XLN(1,  7)  Number  of  contacts; 

XLN(I,  8)  Sum  of  squares  of  down  times; 

XLN(I,  9)  Time  at  which  first  detection  began; 

XLN ( I ,  10)  Time  at  which  last  detection  began; 

XLN(I ,  11)  Time  at  which  current  contact  began  (zero  If  currently 

out  of  contact); 


NTR6TS  Number  of  aircraft  and  barrier  Input  In  TARGTS  array; 

:  (all  times  In  hours) 

All  outputs  are  written  on  unit  6. 

The  following  Is  computed  for  each  target: 


Maximum  value  of  down  time; 

Time  at  which  maximum  down  time  ended; 
Minimum  value  of  down  time; 

Time  at  which  the  minimum  down  time  ended; 
Total  down  time; 

Number  of  contacts; 

Mean  time  between  contacts; 

Standard  deviation  of  down  times; 
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Percent  of  time  covered; 

Average  revisits  per  day; 

Start  time  of  first  contact; 

Start  time  of  last  contact; 

In  addition  the  minimum,  average,  and  maximum  values  across  all 
targets  are  compiled  for: 

Maximum  time  between  contacts; 

Minimum  time  between  contacts; 

Average  time  between  contacts; 

Standard  deviaton  of  down  time; 

Average  duration  of  contacts; 

Percent  coverage; 

Average  revisits  per  day; 

Average  revisit  Interval; 

METHOD: 

Method  is  implied  by  the  outputs. 
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3.5  UTILITIES 


3.5.1  SUBROUTINE  ANGLED 
PURPOSE: 

Computes  the  angular  separation  between  two  points  In  degrees 

INPUTS: 

A(3),  B(3)  —  Unit  vectors. 


OUTPUT: 

ANG  —  Angular  separation  between  A  and  B  In  degrees. 


METHOD: 

Standard  trigonometry. 
REFERENCES: 

No  external  references 
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3.5.2  SUBROUTINE  BRNG 

PURPOSE:  Given  two  pts.,  A  and  B,  find  the  gearing  from  pt.  A  to  pt.  B 
In  radians  and  In  degrees. 

INPUT  ARGUMENTS: 

ALAT  Latti tude  of  pt.  A  (degrees) 

ALON  Longitude  of  pt.  A  (degrees) 

BLAT  Latti tude  of  pt.  B  (degrees) 

BLON  Longitude  of  pt.  B  (degrees) 

OUTPUT  ARGUMENTS: 

RBRNG  Bearing  from  pt.  A  to  pt.  B  In  radians 

DBRNG  Bearing  from  pt.  A  to  pt.  B  In  degrees 

METHOD: 

Given  two  pts.  A  and  B, 

1)  IF  ALAT  =  90.0  or  BLAT  *  -90°  THEN  the  bearing  Is 
equal  to  PI. 

2)  IF  ALAT  =  -90  of  BLAT  =  90°  THEN  the  bearing  Is 
equal  to  0. 

3)  Otherwise 

a)  Find  the  distance  from  pt.  A  to  pt.  B. 

b)  IF  the  difference  between  ALON  and  BLON  Is 
not  180°,  then  the  bearing  from  pt.  A  to  pt. 
B  Is  defined  as  follows: 

cos (Bearing)  -  sin (BLAT)  -  s1n(ALAT)*cos(DIST) 

cos(ALAY)*s1n(DI$t) 

where 

DIST  Distance  from  A  to  B  in  radians 
Bearing  Bearing  from  pt.  A  to  pt.  B 
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c)  IF  the  difference  between  ALON  and  BLON  equals 
180°  then  the  bearing  Is  defined  as  follows: 


IF  AUT  and  BLAT  are  approximately  equal  then 

If  ALAT  BUT  then  bearing  -  f 
If  AUT  BUT  then  bearing  ■  0 

IF.  AUT  and  BUT  do  not  approximate  each  other  then 

If  AUT  0.  then  bearing  ■  0 
If  AUT<0  then  bearing  ■  1 

REFERENCES: 

Subroutine  BRNG  Is  used  by  COVERAGE  and  PROFILE. 


o 
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3.5.3  SUBROUTINE  CIRSCT 


PURPOSE: 

Given  two  points  on  the  earth,  each  with  a  different  radius,  compute 
the  points  of  intersection. 

INPUTS: 

LATA,  LONA  —  Latitude  and  longitude  on  point  A. 

RADA  —  Radius  of  circle  at  point  A  in  NMI. 

LATB,  LONB  —  Latitude  and  longitude  of  point  B. 

RADB  —  Radius  of  circle  at  point  B  in  NMI. 

OUTPUTS: 

LATl,  L0N1  —  Intersect  point  1. 

LAT2,  L0N2  —  Intersect  point  2. 

I POINT  —  Number  of  intersections  (0  or  2). 


METHOD: 


This  routine  will  never  allow  only  one  point  of  Intersection.  If  the 
two  points  are  located  less  than  0.6  degrees  of  each  other,  then  It  Is 
assumed  there  Is  no  Intersection. 

REFERENCES: 

The  following  subroutines  are  called: 

1.  VEC 

2.  CONES 

3.  ANGLED 

4.  RADECM 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


e 


3-58 


vrw  n 


itMihl  It  1 


*,,***  - 


3.5.4  SUBROUTINE  COMBLK 

PURPOSE: 

To  set  up  constants  and  units  for  COVERAGE. 


3.5.5  SUBROUTINE  CONES 
PURPOSE: 

Computes  the  two  unit  vectors  defined  by  the  Intersection  of  two  given 

cones. 

INPUTS: 


A(3)  —  Unit  vector,  axis  of  cone  1. 

B(3)  —  Unit  vector,  axis  of  cone  2. 

C0S1  —  Cosine  of  the  half  angle  for  cone  1, 
C0S2  —  Cosine  of  the  half  angle  for  cone  2. 


OUTPUTS: 


VI (3),  V2(3)  —  The  two  unit  vectors  defined  by  the  Intersection  of 
the  two  cones.  If  the  two  cones  do  not  Intersect,  then  Vl(l)  Is  set  to  2.0. 

METHOD: 


Standard  trigonometry 
REFERENCES: 


No  external  references 


.5.6 


SUBROUTINE  OIS 


PURPOSE:  Given  two  pts.,  A  and  B,  find  the  distance  from  pt.  A  to  pt.  B 
In  radians ,  degrees,  and  nautical  miles. 

INPUT  ARGUMENTS: 

ALAT  Lattitude  of  pt.  A 

ALON  Longitude  of  pt.  A 

BLAT  Lattitude  of  pt.  B 

I OPT  *  1  Implies  coordinates  of  A  and  B  are  Input  In  degrees 

*  2  Implies  coordinates  of  A  and  B  are  Input  In  radians 

OUTPUTS  ARGUMENTS: 

RDIST  Distance  from  pt.  A  to  pt.  B  In  radians 

DDIST  Distance  from  pt.  A  to  pt.  B  In  degrees 

DISTNM  Distance  from  pt.  A  to  pt.  B  In  nautical  miles 

METHOO: 

If  pt.  A  and  pt.  B  are  two  points  on  a  unit  sphere,  then  the 
arc-length  distance  between  them  Is  defined  as  follows 

cos(DIST)  =  s1n(ALAT)s1n(BLAT)  +  cos(ALAT)cos(BLAT)cos(ALON-BLON) 
where  0°<DIST<180° 

output  DIST  In  radians,  degrees,  and  nautical  miles. 

REFERENCES: 

Subroutine  DIST  Is  used  by  COVERAGE  and  PROFILE. 


3.5.7  SUBROUTINE  FLALOB 


PURPOSE:  Given  the  coordinates  of  a  pt.  A,  the  bearing  from  pt.  a  to 
another  pt.,  B,  find  the  coordinates  of  pt.  B. 

INPUT  ARGUMENTS: 

ALAT  Lattitude  of  pt.  A 

ALON  Longitude  of  pt.  A 

OIST  Distance  from  pt.  A  to  pt.  B 
BRNG  Bearing  from  pt.  A  to  pt.  B 

IOPT  -  1  implies  that  all  Inputs  and  outputs  are  In  degrees 
a  2  Implies  that  all  Inputs  and  outputs  are  In  radians 

OUTPUT  ARGUMENTS: 

BLAT  Lattitude  of  pt.  B 

BLON  Longitude  of  pt.  B 

METHOD: 

The  coordinates  of  pt.  B  can  be  defined  as  follows: 

sln(BLAT)  *  cos(BRNG)cos(ALAT)s1n(DIST)  +  s1n(ALAT)cos(DIST) 

BLON  -  SIGN(DLONG,  BRNG) 

where  DLONG  Is  the  difference  between  ALON  and  BLON 
cos(DLONG)  = 

BLON  Is  then  adjusted  such  that 
0<I BLON | <180° 

REFERENCES: 

FLALOB  Is  used  by  COVERAGE. 
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3.5.8  Subroutine  HISTGM 


PURPOSE: 

HISTGM  is  a  utility  routine  responsible  for  generating  a 
horizontally  arrayed  histogram  bar  chart  on  a  character  printer  (as 
opposed  to  a  graphics  display  device).  Both  self-scale  and  user-scale 
options  are  available.  Both  percent  and  absolute  levels  are  printed. 

INPUT  ARGUMENTS: 

NCOUNT  Total  of  all  frequency  counts; 

I FREQ  Array  of  frequency  counts  (rows) ; 

NROW  Number  of  bins  or  rows; 

MAX  Scaling  term: 

>  0  for  the  maximum  scale  (frequency  for  any  row); 

±  0  for  self-scale; 

LABEL  Array  of  eight  character  per  row  alphanumeric  row  labels; 

I UNIT  The  FORTRAN  logical  unit  number  of  the  file  on  which  the 

histogram  Is  to  be  written. 

OUTPUT  ARGUMENTS: 

None 

METHOD: 

If  MAX  Is  less  than  or  equal  to  zero,  the  routine  finds  the  maximum 
frequency  for  all  the  rows. 

All  frequencies  will  be  scaled  by  50. /maximum. 

The  histogram  Is  a  bar  chart  whose  length  Is  computed  by: 

J  -  TFIX(FLOAT(IFREQ( I ))*(50. /Maximum)  +  .5) 

J  *  MIN0(50,  MAX0(0fJ)) 


3.5.9 


SUBROUTINE  JDORGF 


PURPOSE: 

To  convert  a  date  from  Julian  to  year,  month,  day,  hours,  minutes, 
seconds  or  from  year,  month,  day,  hours,  minutes,  seconds  to  julian  date. 


3.5.10  SUBROUTINE  LINTRP 


PURPOSE: 

LINTRP  is  a  utility  routine  responsible  for  performing  a  bounded 
linear  interpolation  of  the  value  of  a  dependent  variable  given  the 
Independent  variable  and  a  table  of  dependent  and  Independent  variables. 
The  Independent  variable  must  be  monltonically  Increasing. 

INPUT  ARGUMENTS: 


XI  Independent  variable 

X  Table  of  Independent  variables 

Y  Table  of  dependent  variables 

N  Length  of  tables 

OUTPUT  ARGUMENTS: 


Y1  Dependent  variable  associated  with  XI. 

METHOD: 


Yo  If  XI  i  Xo 

Yl  -  Y1  +(Y1  -  Y1-1  (XI  -  X1-1)/(X1  -  X1-1)  for  X1-1  <  XI  <  XI 
Yn  if  XI  >  Xn 
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3.5.11 


SUBROUTINE  MAGN 


PURPOSE: 

To  compute  the  magnitude  of  vector  A. 

INPUTS:  A  A  three  element  array. 

OUTPUTS: 

B  Magnitude  of  A; 

C  Magnitude  of  A  squared; 

METHOD: 

Magnitude  of  vector  A  Is  the  square  root  of  the  sum  of  the  squares  of 
the  three  elements  of  A. 

REFERENCES: 


MAGN  Is  used  by  COVERAGE 
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3.5.12  SUBROUTINE  RADECM 
PURPOSE: 

Computes  the  right  ascension  and  declination  of  a  given  vector,  and 
Its  magnitude. 

INPUTS: 


ALPHA  —  Right  Ascension  In  degrees.  • 
DELTA  —  Declination  In  degrees. 

VMAG  —  Magnitude  of  V. 

METHOD: 

Standard  trigonometry 
REFERENCES: 

No  external  references. 
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3.5.13 


SUBROUTINE  RMMUL 


PURPOSE: 

To  find  the  product  of  two  matrices. 

INPUTS: 

A 
B 

OUTPUTS: 

C  Product  of  matrices  A  and  B; 

(Matrix  of  size  L  x  N) 

REFERENCES: 

RMMUL  Is  used  by  COVERAGE. 


Matrix  of  size  L  x  M; 
Matrix  of  size  M  x  N; 


3.5.14  SUBROUTINE  SPHTRN 

PURPOSE:  Given  the  coordinates  of  a  point  P  in  the  old  system  and  the 
coordinates  of  a  new  pole  in  the  old  coordinate  system,  find 
the  coordinates  of  point  P  in  the  new  system. 

INPUT  ARGUMENTS: 

PHI  Longitude  of  pt.  in  old  system 
TH  Colatitude  of  pt.  In  old  system 

PHIO  Longitude  of  the  new  pole 
THO  Colatitude  of  the  new  pole 
PPHIO  Bearing  of  the  azimuthal  reference 

OUTPUT  ARGUMENTS: 

PPHI  Longitude  of  pt.  In  the  new  system 
PTH  Colatitude  of  pt.  in  the  new  system 

METHOD: 

The  new  azimuthal  reference  Is  0°  0.  Therefore  the  colatitude 
of  the  pt.  P  In  the  new  system  can  be  defined  as  follows: 

cos(PTH)  =  cos (THO)cos (THO)cos (TH)+s 1 n (THO) sin (TH)cos (PHIO) 

where  PTH  Is  the  colatitude  of  pt.  P  In  the  new  system. 

sin(PPHI)  «  -s1n(PHI0)s1n(TH)/s1n(PTH) 

REFERENCES: 

SPHTRN  Is  used  by  COVERAGE. 
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3.5.15  Subroutine  STHIST 


PURPOSE: 

Tills  utility  routine  Is  responsible  for  computing  the  mean  and 
standard  deviation  of  a  set  of  data  and  outputting  the  minimum  and 
maximum,  total  count,  mean  and  standard  deviation  and  frequency  histogram. 

INPUT  ARGUMENTS: 

NCOUNT  Total  of  all  frequency  counts,  N; 

I FREQ  Array  of  frequency  counts; 

NROW  Number  of  rows  or  bins  Into  which  frequency  has  been 
divided; 

MAX  Explicit  scale  for  histogram  If  greater  than  0; 

LABEL  Eight  character  row  labels  for  histogram; 

I UNIT  FORTRAN  logical  unit  number  for  output; 

XMIN  Minimum  value  of  X; 

XMAX  Maximum  value  of  X; 

SUM  SUM  of  X; 

SUM2  SUM  of  X  squared; 

OUTPUT  ARGUMENTS: 

SUM  Mean  value  of  X 

SUM2  Standard  deviation  of  X 

METHOD: 

The  mean  and  standard  deviation  of  X  are  computed  as  follows  for 
NCOUNT  great  than  1: 


SUM  -  SUM/FL0AT( NCOUNT) 

SUM2  -  SQRT((SUM2  -  FLOAT(NCOUNT)*SUM*SUM)/FLOAT(NCOUNT-1) ) 


it  *  \** 

,1  V;> 


The  routine  writes  the  NCOUNT,  minimum,  maximum,  mean  and  standard 
deviation  on  device  IUNIT.  It  then  uses  routine  HISTGM  to  generate  a 
frequency  histogram. 

REFERENCES: 

STHIST  Is  called  by  STATPT  routine  and  It  calls  HISTGM. 


3.5.16 


SUBROUTINE  SUN 


PURPOSE: 


To  compute  the  orbital  elements  of  the  sun  for  a  given  Julian  date. 

INPUTS: 

TOTJD  JULIAN  DATE 

OUTPUTS: 


0ES(6)  Orbital  elements  of  the  sun; 

REFERENCES: 

SUN  Is  used  by  coverage. 


G 


BaraBWMaMUffiJfja 


b.’ v'  •_*.■«%*•.  *.*  »• 


3.5.17  SUBROUTINE  TRI 


PURPOSE: 

61ven  either  two  sides  and  an  angle  of  a  triangle  or  two  angles  and  a 
side  of  a  triangle,  an  option  flag,  and  an  ambiguity  flag,  find  the  remaining 
sides  and  angles  of  the  triangle. 


INPUT  -  OUTPUT  MATRIX: 


DEFINE  P,  Q,  R  -  angles  (degrees) 

X,  Y,  Z  -  sides  (degrees)  such  that 
x  Is  opposite  P, 
y  Is  opposite  Q,  and 
z  Is  opposite  R. 

L  -  option  Flag 

M  -  ambiguity  Flag 


To  OUTPUT: 

ANGLES,  SIDES  (RAD) 

R(QUAD  1  or  2),  X,  Y  1 

R(QUAD  3  or  4),  X,  Y  1 

Y(QUAD  2  or  3),  R,  Z  2 

X(QUAD  1  or  2),  0,  R  3 

X(QUAD  3  or  4),  Q,  R  3 

Q(QUAD  1  or  4),  R,  Z  4 

Q(QUAD  2  or  3),  R,  Z  4 

P(QUAD  1  or  2),  Q,  R  5 

P(QUAD  3  or  4),  Q,  R  5 

X(QUAD  1  or  2),  Y,  Z  6 

X(QUAD  3  or  4),  Y,  Z  6 


ANGLES,  SIDES  (RAD) 

1  P,  Q*  Z 

2  P,  Q,  Z 

2  P,  Q,  X 

1  P,  Y,  Z 

2  P,  Y,  Z 

1  P,  X,  Y 

2  P,  X,  Y 

1  X,  Y,  Z 

2  X,  Y,  Z 

1  P.  Q,  R 

2  P,  Q,  R 


INPUT: 

L  M 


METHOD: 


Simple  trigonometry 


3.5.18  SUBROUTINE  VEC 


i 


PURPOSE: 

Computes  the  unit  vector  corresponding  to  a  given  right  ascension  and 
declination  In  degrees. 

INPUTS: 

ALF  —  Right  ascension  In  degrees. 

DEL  —  Declination  In  degrees. 

OUTPUTS: 

V(3)  —  Unit  vector  corresponding  to  ALF  and  DEL. 

METHOD: 

Standard  trigonometry 
REFERENCES: 

No  external  references. 


3.5.19 


SUBROUTINE  ZERRV 


3.5.20  SUBROUTINE  ZERIV 


PURPOSE: 

To  Initialize  an  Integer  array. 


INPUTS: 


An  Integer  array  of  dimension  N; 
Dimension  of  IA; 


OUTPUTS: 


IA(I)  *  0  where  11  I  <  N; 


REFERENCES: 

This  routine  Is  used  by  COVERAGE  and  AAWSIM. 
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lows  are  utilized  In  the  program 


the  cosine  of  an  angle,  find  the  angle  In 


